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1 .  INTRODUCTION 

The  problem  of  electromagnetic  coupling  to  a  cavity 
through  an  aperture  is  an  old  problem.  For  an  extensive  list 
of  references  see  Cl  3.  A  number  of  problems  are  solved  by 
the  so-called  generalized  network  formulation  for  aperture 
problems  C23— C43.  In  this  method;  the  aperture  is  closed  by 
a  perfect  conductor;  a  sheet  of  unknown  magnetic  current  M 
is  placed  on  one  side  of  the  shorted  aperture  and  -M  is 
placed  on  the  other  side.  This  insures  continuity  of  the 
tangential  electric  field  across  the  aperture.  Next;  the 
tangential  magnetic  field  is  forced  to  be  continuous  across 
the  aperture;  resulting  in  an  integral  equation  for  M.  M  is 
found  numerically  by  following  the  method  of  moments  C5D . 
Unfortunately;  this  method  can  be  easily  implemented  only 
when  the  aperture  is  in  a  perfectly  conducting  infinite 
plane  so  that  image  theory  C6;Sec.3-43  may  be  applied  to 
find  the  fields  due  to  M  C33;C43. 

When  the  aperture  is  in  a  finite  curved  surface;  the 
fields  due  to  M  (or  -M)  radiating  in  the  presence  of  the 
complete  conducting  surface  are  difficult  to  find.  The 
complete  conducting  surface  is  the  conducting  surface  with 
its  aperture  shorted.  The  present  report  gives  an  accurate 
solution  to  such  problems.  We  recognize  that  the  problem  of 
obtaining  the  electromagnet ic  field  due  to  M  is  a  scattering 
problem  in  which  the  impressed  source  is  on  the  scatterer.  A 


similar  static  problem  was  solved  in  C71  by  a  method  in 
which  a  pseudo-image  was  used  in  order  to  improve  the 
numerical  accuracy.  For  reference,  we  call  this  method  the 
pseudo-image  method.  With  this  in  mind,  we  approach  the 
problem  by  following  the  generalized  network  formulation 
method  until  the  point  of  finding  the  fields  due  to  hi 
radiating  in  the  presence  of  a  complete  conducting  surface. 
We  then  solve  this  scattering  problem  by  placing  an  electric 
current  J  on  the  complete  conducting  surface,  introducing 
the  pseudo-image  as  in  C71,  writing  the  electric  field 
integral  equation  for  J  and  finally  solving  this  equation 
for  J  by  the  method  of  moments  C53.  As  shown  in  the 
following  sections,  the  theory  developed  in  this  way  works 
very  well  for  a  non-resonant  cavity  formed  by  the  complete 
conducting  surface.  However,  it  fails  when  the  cavity  is 
resonant.  Further  study  is  needed  to  treat  this  special 
case. 

Section  2  states  and  formulates  the  general  problem 
shown  in  Fig.l.  The  theory  developed  applies  to  both  two  and 
three  dimensional  problems.  In  a  two  dimensional  problem, 
all  quantities  are  invariant  in  one  direction,  for  instance, 
the  z-direction.  Section  3,  which  is  essentially  the 
application  of  the  theory  developed  in  Section  3,  solves  the 
problem  of  electromagnetic  scattering  from  an  infinitely 
long,  infinitesimally  thin,  perfectly  conducting  cylindrical 
surface  with  an  infinitely  long  slot  illuminated  by  a 
uniform  TM  plane  wave.  Section  4  states  the  numerical 


results  obtained  for  the  problem  of  Section  3.  Various 
methods  are  compared  to  show  how  well  the  pseudo— image 


method  works.  Appendix  A  evaluates  some  of  the  complicated 
integrals  appearing  in  Section  3.  Appendix  B  summarizes  a 
method  for  the  aperture  problem  in  which  the  aperture  is  not 


closed  by  a  perfect  conductor  and  only  an  unknown  electric 
current  J  is  placed  on  the  conducting  surface.  This  method 
fails  to  give  accurate  fields  inside  the  cavity  when  the 
aperture  becomes  very  small. 
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a .  FORMULATION  OF  THE  PROBLEM 

a . 1  Statement  of  the  Problem 

The  problem  to  be  considered  is  shown  in  Fig.l  where 
electromagnetic  coupling  occurs  through  an  aperture  of 
arbitrary  shape  in  a  perfectly  conducting  body  which  may 
have  finite  thickness.  If  the  aperture  were  closed,  there 
would  be  a  cavity.  Some  additional  complete  perfectly 
conducting  bodies  may  exist  either  inside  or  outside  the 
cavity,  or  both.  To  keep  the  formulation  of  the  problem 
simple,  we  consider  the  case  where  only  one  additional 
complete  conducting  body  is  located  outside  the  cavity. 
Other  cases  are  just  straightforward  extensions  of  the  case 
under  consideration.  A  homogeneous  medium  with  permeability 
H  and  permittivity  €  fills  all  space  outside  the  perfectly 
conducting  bodies.  The  perfectly  conducting  body  with  an 
aperture  is  specified  by  its  internal  surface,  denoted  by 
S* ,  and  its  external  surface,  denoted  by  Sm.  (Of  course,  if 
the  body  has  zero  thickness,  then  the  internal  and  external 
surface  will  be  indistinguishable.)  The  additional 


perfectly  conducting  body  is  specified  by  the  surface  S« .  It 
does  not  matter  whether  this  body  is  hollow  or  solid 
because,  due  to  the  fact  that  Sm  is  closed,  no  field  can 
penetrate  it.  The  cross  sections  of  the  conducting  bodies 
are  shown  in  Fig.l.  The  problem  is  two  dimensional  if 
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everything  is  invariant  in  the  z— direction.  In  this  case  the 
conducting  bodies  are  infinitely  long  cylinders  and  the 
aperture  is  an  infinitely  long  slot.  We  formulate  the 
problem  without  mentioning  whether  it  is  a  three-  or  two- 
dimensional  problem. 

Next *  we  define  the  symbols  appearing  in  Fig.l. 
( J M*"***  >  denotes  the  known  impressed  electric  current 
source  and  impressed  magnetic  current  source.  (Elr,e  *Hlr,c> 
denotes  the  fields  that  would  exist  if  ( )  were  to 
radiate  into  unbounded*  homogeneous  space  with  ()»*€) 
everywhere.  Region  b  denotes  the  space  inside  the  cavity* 
and  region  a  denotes  the  space  external  to  the  closed 
surface  that  consists  of  S.  and  the  aperture. 
(Ekr^+E® denotes  the  total  field  in  region  a*  and 
( E*3  *  HP )  denotes  the  total  field  in  region  b.  The  objective 
is  to  determine  these  total  fields.  In  turn*  related 
quantities  such  as  power  gain  and  polarizability  can  be 
easily  determined.  Also*  we  want  to  investigate  how  Bethe 
hole  theoryCSD » C9D  should  be  modified  when  a  small  aperture 
(compared  to  wavelength)  is  in  a  finite  perfectly  conducting 
cavity  instead  of  in  an  infinite  ground  plane.  Finally,  we 
state  the  boundary  conditions  required  for  solving  this 
problem.  First,  both  the  tangential  electric  field  and  the 
tangential  magnetic  field  must  be  continuous  across  the 
aperture.  Second*  the  tangential  electric  field  on  the 
perfectly  conducting  surfaces*  i.e.,  on  Sv*  S*  and  S«  must 
be  zero.  With  this  in  mind*  we  proceed  to  the  next  section. 


‘>1 


w 


m 

a 

$•] 


i 


»w 


aperture 


region  b 


(**) 


£b.Hb 


region  a 

(/i,€  ) 


«S 


Original  problem:  Electromagnetic  coupling  through  an 
arbitrary  shaped  aperture  in  a  conducting  body.  (Only 
a  typical  cross  section  is  shown;  the  shaded  area  denotes 
a  perfect  conductor.)  Impressed  sources  and 

radiate  in  the  presence  of  surfaces  S^. ,  S^,  and  S  . 
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S . 3  Derivation  of  the  Operator  Equation 

The  situation  considered  in  Fig.l  is  rather 
complicated.  The  approach  that  we  take  is  to  decompose  the 
problem  into  two  parts.  The  point  of  departure  is  the 
equivalence  principle  £6*Sec.3-51.  The  aperture  in  Fig.l  is 
first  closed  by  a  perfectly  conducting  flat  plate.  (If,  in  a 
three  dimensional  case*  it  is  not  possible  to  close  the 
aperture  by  a  flat  plate*  then  special  treatment  is  needed.) 
Then  the  magnetic  current  sheet  M  is  placed  on  the  left-hand 
side  of  the  flat  plate  and  -M  is  placed  on  the  right-hand 
side  of  the  flat  plate.  The  M  is  defined  as 

rt=  E  x  n  (3-1) 

where  n  is  the  unit  normal  vector  on  the  aperture*  pointing 
toward  region  a*  and  E  is  the  electric  field  in  the  aperture 
in  Fig.l.  In  this  way*  the  boundary  condition  that  the 
tangential  electric  field  is  continuous  across  the  aperture 
is  satisfied.  The  combination  of  S«  and  the  flat  plate  is 
called  S.  ;  and  the  combination  of  Si  and  the  flat  plate  is 
called  .  (Superscript  sc  stands  for  short  circuit  since 
the  closing  of  the  aperture  can  be  viewed  as  the  short 
circuiting  of  the  aperture).  The  equivalent  problem  is  shown 
in  Fig. 3.  Before  the  problem  in  Fig. 3  is  actually  equivalent 
to  the  problem  in  Fig.l*  we  have  to  impose  continuity  of  the 
tangential  magnetic  field  across  the  aperture*  as  required 
in  the  original  problem.  Namely* 
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H«(J  **"*»,  M“"*>  >+tU<0jM>=H«  (0,-M) 


on  A 


( 2—2  > 


where  A  denotes  the  flat  plate  shorting  the  aperture*  The 
subscript  t  in  (2-2)  denotes  the  component  tangent  to  A.  In 
(2-2),  H*(  -  ,  - )  is  the  magnetic  field  operator  in  region  a 
with  the  sources  indicated  in  the  parentheses  (-,-) 

me 

radiating  in  the  presence  of  S.  and  S«.  H*»(-,->  is  the 
magnetic  field  operator  in  region  b  with  the  sources 
indicated  in  the  parentheses  (-,-)  radiating  in  the  presence 

me 

of  S*.  Actually, the  '‘a”  fields  on  the  left-hand  side  of  (2- 
2)  are  evaluated  not  exactly  on  A  but  immediately  to  the 
left  of  any  electric  current  that  could  flow  on  the  flat 
plate.  Moreover,  the  "b"  field  on  the  right-hand  side  of  (2- 
2)  is  evaluated  not  exactly  on  A  but  immediately  to  the 
right  of  any  electric  current  that  could  flow  on  the  flat 
plate.  Each  magnetic  field  in  (2-2)  has  two  arguments.  The 
first  argument  is  an  electric  current  source  and  the  second 
argument  is  a  magnetic  current  source.  Eguation  (2—2)  states 
that  the  tangential  magnetic  field  must  be  continuous  across 
the  aperture. 

The  impressed  sources  Jlm**  and  M*mF’  and  the  boundary 
conditions  on  the  fields  on  A,  S*.  ,  S.,  and  S«  in  Fig. 2  are 
the  same  as  in  Fig.l.  It  is  evident  from  C 10, Theorem  II  that 
the  electric  field  in  an  infinite  region  outside  a  closed 
surface  is  uniquely  specified  by  its  impressed  sources  in 
the  region  and  its  tangential  components  on  the  surface. 
Consequently,  the  fields  in  Fig. 2  are  identical  to  those  in 


Fig.l.  Comparing  Fig.l  with  Fig. 2,  we  find  that 

<E‘"‘=+E-,H1"«=+t?-)=  (2-3) 

CE-(  Jk™»*»Mfc",r*)+i*(0»M)  *H-(  J*"*»  ,  M 1  m(*  >  +H*  ( 0  >  M )  > 

in  region  a 

)  =  >  in  region  b  (2-4) 

The  electric  fields  on  the  right-hand  sides  of  (2-3)  and  (2— 
4)  correspond  to  the  magnetic  fields  in  (2-2).  E* (  -  ,  -  )  is 
the  electric  field  operator  in  region  a  with  the  sources 
indicated  in  the  parentheses  radiating  in  the  presence  of 

me 

S„  and  S*.  E* ( - f - >  is  the  electric  field  operator  in  region 
b  with  the  sources  indicated  in  the  parentheses  radiating  in 
the  presence  of  S».  .  For  example,  Eto(0,-M)  is  the  electric 
field  in  region  b  due  to  -M  placed  on  the  right-hand  side  of 
the  flat  plate  that  shorts  A. 

To  find  the  unknown  magnetic  current  M,  we  appeal  to 
(2—2).  Rearranging  the  terms  in  (2-2)  and  using  the 
linearity  of  the  operators,  we  obtain 

m  to  m 

-Ht,(0,M)-H*  (0,M)=H.%<  Jlm**  ,Mkm**)  on  A  (2-5) 

2 . 3  Numerical  Solution  of  the  Operator  Equation 

To  find  M  using  (2-5),  we  first  approximate  M  by 


where  M 


are  known 


(to  be  chosen)  vector  expansion 


functions)  and  Vn  are  scalar  coefficients.  Substituting  (2- 
6)  into  (2-5) »  we  have*  by  the  linearity  of  the  operators  in 
(2-5)  , 

«  b  m 

Z  v„(— H*(0,M„)— H*(0,M„)  )=H%(  (2-7) 


Next*  we  define  a  symmetric  product  as 


A,B  >=  J 


A  •  B  ds 


(2-8) 


where  A  and  B  are  vector  functions*  S  is  the  surface  on 
which  they  are  defined*  and  ds  is  the  differential  element 
of  surface  area  on  S.  Taking  the  symmetric  product  of  (2-7) 
with  rim,  m=l,2,-*-,  we  obtain 


CY-]-»-CY*»l 


V  =  I*™ 


(2-9) 


where  V  and  I*’’®  are  column  vectors.  The  n*H  element  of  V  is 

■4  %  nc 

The  m**  element  of  I1"®  is  Im  given  by 


Im  =<Mm,H*(J‘-*»,M*"*»)> 

Moreover  *  CY*3  and  CY**]  are  square  matrices  whose  mn*h 
elements  are  given  by 


(2-10) 


-I, 

-L 


Mm  •  H%(0,M„)  ds 


Mm  *  H%(0,Mr»)  ds 


(2-11 ) 


(2-12) 


12 


s- 


In  <2-11)  H~(0»Mn)  is  evaluated  on  the  region  a  side  of  A. 

Similarly,  in  (2-12),  Hto(0,M„)  is  evaluated  on  the  region  b 
side  of  A.  If  CY*],  CY**D,  and  I **■*«=  are  found,  V  can  be  found 
by  solving  the  set  of  linear  equations  given  in  (2-9).  A 
standard  computer  routine  for  this  purpose  is  available. 

Substituting  (2-6)  into  (2-3)  and  (2-4),  we  relate  each 
field  component  in  those  two  equations  to  the  fields  due  to 
the  known  expansion  functions  (sources).  Now  the  problem  of 
finding  the  fields  due  to  the  unknown  source  M  is  reduced  to 
the  problem  of  finding  the  fields  due  to  known  sources 
radiating  in  the  presence  of  complete  perfect  conductors. 
The  fields  in  each  region  can  be  written  as  follows:  In 
region  a. 


E*™=+E«  =  E-<  )+Z  V„  E»(0,Mr.) 


Hir.e+y.  =  H—  >+Z  V„  H-<0,M„) 


In  region  b; 


(2-13) 


(2-14) 


E*»  =  -Z  V„  Eb(0,f«U) 


HP*  =  - Z  Vn  H*>  (0,M„> 


(2-15) 


(2-16) 


The  fields  on  the  right-hand  sides  of  equations  (2-13)-(2- 
16)  have  the  same  meaning  as  in  (2-3)  and  (2-4).  The  fields 
on  the  left-hand  sides  of  <2-13)-(2-16)  are  those  in  the 
original  problem  Fig.l. 

In  the  following  sections,  we  evaluate  the  fields  on 
the  right-hand  sides  of  (2-13)-(2-16)  one  by  one.  The 


13 


results  obtained  are  not  only  substituted  into  (2-13)-(2-16) 
but  also  into  (2-10>— (E-1S>  to  find  the  coefficients  of  the 
magnetic  current  expansion  functions*  i.e.*  V„,  n=l,2,-»«. 

2.4  Short-Circuit  Field  due  to  the  Impressed  Sources 

(E-(  ,H-(J4m*»,M4m*»)  )  in  (2-13)  and  (2-14)  are 

the  fields  in  region  a  due  to  the  impressed  sources 
<  J4*"**,  M4*"** )  radiating  in  the  presence  of  S_  and  S«  (See 
Fig. 2).  These  fields  can  be  viewed  as  the  sum  of  the  fields 
due  to  the  impressed  sources  radiating  in  an  unbounded 
homogeneous  space  with  (}*»€)  and  the  fields  due  to  the 

mcz 

electric  current  induced  on  S»  and  S_  radiating  also  in  the 
unbounded  homogeneous  space  with  (H>€).  Namely*  in  region  a* 


E*(  »M *«•*»)  =  E4”=+E(J*’*,0) 

H*(Jl"">,M‘”»)  =  Hlr~r+H<  J-M  ,0) 


(2-17) 

(2-18) 


where  (E4r,<s * Hine )  is  defined  in  the  last  paragraph  of 
section  (2.1)*  (E(  J***,  0)  ,H<  Jmm  *0)  )  are  the  fields  due  to 
induced  current  JmM  on  Sm  and  S«,  radiating  in  an  unbounded 
homogeneous  space  with  (p*€). 

me 

□n  the  surfaces  of  the  perfectly  conducting  bodies  S«, 
and  S« 

•  1  ne 

I*  <  J *"•**, =  E*  +E*(J— ,0>  =  0  (2-19) 

where  the  subscript  t  denotes  the  component  tangent  to 

•c 

either  S.  or  S«.  Rearranging  (2-19),  we  obtain 
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tV 


V 

IL  ■  » 


w 


— E%(  J***  »0)  ■  E* 

«C 

on  S»  and  S« 

The  electric  current  J***  is  approximated  by 


J— ■=  2  Ij  J3 

j 


(2-20) 


(2-21 ) 


where  <Jj>  are  expansion  functions  that  are  chosen  later 
when  particular  impressed  fields  are  given.  The  function  J,» 
is  on  either  S.  or  S«  and  is  tangent  to  whichever  surface  it 
is  on.  Seeking  to  determine  the  unknown  coefficients  Ij  ,  we 
substitute  (2—21)  into  (2—20)  and  take  the  symmetric  product 
of  (2—20)  with  Jii  i=l,2,---,  to  obtain 


C2*3  I-**  = 


(2-22) 


where  I*M  and  V-**  are  column  vectors.  The  j**'  element  of  I"** 

-*  *  N 

is  Ij  .  The  i*r'  element  of  V***  is  V»  given  by 


-  <  3%  , E 1  > 


(2-23) 


Furthermore^  CZ*3  is  a  square  matrix  whose  i  j element  is 
given  by 


=  <J»  *— E(  J 


I,0)>  =  -|  j! 


J*  '  *0)  ds 

sl%s« 


( 2-24 ) 


where  |  denotes  that  the  surface  integral  is  performed  on 
S.  +S. 

»«=  m, 

both  S.  and  Sa.  Note  that  Elr,c:  and  E(Jj,0)  are  due  to  their 
corresponding  sources  radiating  in  empty  space.  They  can  be 
found  by  the  well  known  formulas  in  empty  space.  In  turn  (2- 


mmmmm 
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23)  and  <2-24)  can  be  computed.  With  E2*D  and  V***  available, 

MM 

<Ij  >  can  be  found  by  solving  the  matrix  equation  (2—22). 
Substituting  (2-21)  into  (2-17)  and  (2-18),  we  obtain,  in 
region  a, 


E-  <  J  ,  M*«e  )  =E  *™=  +  2  Ij  E(Jj,0)  (2-25) 

a 

+  Z  Ij  H  ( J  j  ,  0 )  (2-26) 

As  introduced  in  Fig. 2,  the  fields  on  the  left-hand  sides  of 
(2—25)  and  (2-26)  exist  only  in  region  a  so  that  (2-25)  and 
(2-26)  are  only  valid  in  region  a  even  though  (E1"*5  ,H*-™=) 

m  m 

and  (E ( J j , 0) ,H( Jj ,0) )  are  valid  everywhere  in  space. 


2-5  Short-Circuit  Fields  due  to  the  Magnetic  Current 

Expansion  Functions 

a)  Finding  <E«< 0,M„ ) ,H-< 0, M„ ) ) 

The  field  (E*(  0,  M„  >  ,H*<  0,  Mr,)  )  appearing  in  (2-13)  and 
(2—14)  is  the  field  in  region  a  due  to  the  magnetic  current 
sheet  Mr.  placed  on  the  region  a  side  of  the  shorted  aperture 
and  radiating  in  the  presence  of  S_  and  S. .  Since  S.  is 
perfectly  conducting,  the  region  inside  it  is  completely 
isolated  from  region  a.  Therefore,  any  medium  and  any  source 

me 

can  be  placed  inside  S»  without  affecting  the  field  in 

•C 

region  a.  We  fill  the  region  inside  S.  with  the  medium 
characterized  by  <H>€)  and  place  the  magnetic  current  sheet 
M„  on  the  right-hand  side  of  the  shorted  aperture  as  shown 
in  Fig. 3.  The  magnetic  current  sheet  on  the  right-hand  side 


of  the  shorted  aperture  in  Fig. 3  is  called  the  pseudo-image 
of  the  one  on  the  left-hand  side  of  the  shorted  aperture. 
Now ,  the  electromagnetic  field  (E*(0,M„)  >H*(0,M„)  )  in  region 
a  is  written  as 


E*<0»M«)=E(0»2Mr,>+E(O„,0)  (3-27) 

H“(0,M„)=H(0,2Mr,)+H(2„,0)  (2-28) 


The  fields  on  the  right-hand  sides  of  (2-27)  and  (2-28)  are 
due  to  their  correspond ing  sources  given  in  the  associated 
parentheses  radiating  in  the  unbounded  homogeneous  medium  as 

m 

on  the  right-hand  sides  of  (2-17)  and  (2-18).  Here*  V.-,  is 
the  electric  current  induced  on  S„  and  S.  by  2M„  where  2M„ 
represents  the  combination  of  the  two  magnetic  current 
sheets  in  Fig. 3.  Now,  adjusts  itself  so  that 

m 

E*(0,lin)  =0  (2-29) 

on  Se  and  Sa 

where  subscript  t  denotes  tangential  component.  Combining 
(2-27)  and  (2-29),  we  obtain 
a 

-  E%(3n,0)  =  I*(0,2M„>  (2-30) 

2r»  is  now  approximated  by 

2„  =  E  I 3„  Jj  (2-31 ) 

J 

m 

where  I j  =  are  coefficients  to  be  determined  *  and 

m 

Jj  are  the  same  as  in  (2-21).  Substituting  (2-31)  into  (2- 
30)  and  taking  the  symmetric  product  of  (2-30)  with  each 


Jt,  we  obtain 


Here,  CZ*1  is  given  in  (2-24),  and  I„  and  V„  are  column 

m 

vectors.  The  j**'  element  of  I„  is  I and  the  i **  element 

-*m  m 

of  V„  is  given  by 

=  <3x  ,E*(0,2M„>>  (2-33) 

Because  E(0,2Mr.)  is  the  field  due  to  2Mr,  radiating  in 
unbounded  homogeneous  space,  it  can  be  evaluated  easily. 
Thus,  (2-33)  can  be  computed.  Therefore  In  can  be  found  by 
solving  the  linear  equations  in  (2-32),  i.e.,  is 

completely  determined.  Substituting  (2—31)  into  (2—27)  and 
(2-28),  we  obtain,  in  region  a. 


m  m 

E«(0,M„)=E(0,2Mt,)  +  Z  E(Jj,0)  (2-34) 

H*  (0  ,M„)  =H(  0, 2Mr, )  +  z  I3n  H(J,,0)  (2-35) 

In  (2-33)— (2-35) ,  2M„  represents  the  combination  of  the 
magnetic  current  sheets  placed  on  both  sides  of  the  flat 
plate  that  shorts  the  aperture  in  Fig. 3.  At  any  point  not  on 
the  shorting  plate  E(0,2Mr,)  and  H(0,2M„)  are  simply  twice 
the  fields  due  to  M„  because  the  two  M„’s  are 
infinitesimally  close  to  each  other.  On  the  shorting  plate 
the  normal  electric  field  E„(0,2M„)  and  the  tangential 
magnetic  field  Ht^giyEtl,,)  are  simply  twice  those  due  to  M„, 
but  E%(0,2Mr,)  and  H„(0,2M„)  vanish.  Because  E«(0,2f1„) 
vanishes  on  the  shorting  plate,  the  magnitude  of  the 
electric  current  on  it  is  considerably  reduced.  Hopefully, 


this  reduction  will  improve  the  behavior  of  the  electric 


current  on  the  shorting  plate  and  consequently  a  more 
accurate  numerical  solution  to  the  problem  can  be  obtained. 

b)  Finding  (E*  (0,M„>  ?HP»  (0,M,?>  ) 

<£to<0»Mr?> ?Hto(0?M„) )  in  (£-15)  and  (£-16)  are  the  fields 
due  to  M,-,  placed  on  the  right-hand  side  of  the  shorted 

■c  «c 

aperture  radiating  in  the  presence  of  S*  .  Since  S»  forms  a 
cavity?  will  produce  zero  fields  in  region  a.  With  the 
same  idea  as  in  Section  £.5a?  we  write  for  region  b 


ib<0»M«)asE(0,£M„)+E<2^»0) 

to 

hf  <0?Mr,)=H(0,£Mri)+H(2r?»0) 


( £-36 ) 
( £-37 ) 


where  <E<0»2Mr.)  ,H(0»£M„)  )  is  the  electromagnetic  field  due 
to  air,  radiating  in  unbounded  homogeneous  space  where  £tir> 
represents  the  combination  of  the  two  magnetic  current 
sheets  shown  in  Fig. 4.  There?  the  placed  on  the  right- 

hand  side  of  the  shorted  aperture  is  the  original  M„,and  the 
Mn  placed  on  the  left-hand  side  of  the  shorted  aperture  is 

to 

called  the  pseudo  — image  of  Mr?.  Dr?  is  the  electric  current 

•e  b  b 

induced  by  on  Si  .  (  E(D„,  0)  ?  H(Dr??  0)  )  is  the 

to 

electromagnetic  field  due  to  Dr?  radiating  in  homogeneous 

to 

space.  Now?  Dr?  adjusts  itself  such  that 


1.(0, Mr?)  =  0 


on  S  t 


( £-38 ) 


where  subscript  t  denotes  tangential  component.  Combining 
(£-36)  and  (£-38)?  we  obtain 


■•WWW! 


(QMn) 


shorted  ^  l  re»onb 
aperture  A  |\(  /x,*) 


The  xy  cross  section  of  a  situation  in  which  (E  (0,  M  ), 
b  —  -n 

H  (0,  M  ))  can  exist.  The  magnetic  current  sheet  immedi 

ately  to  the  left  of  the  shorted  aperture  is  called  the 
pseudo-image  of  the  one  immediately  to  the  right. 


?l» 


-E*<5r,,0>=Eti  (0?2M„)  on  S,. 


( H-39 ) 


is  next  approximated  by 


3n  =  2  Iln  Jj 


( 2-40 ) 


where  Jj  is  an  expansion  function  on  S*  .  Substituting  (2- 
40)  into  (2—39)  and  taking  the  symmetric  product  of  (2—39) 

fc> 

with  Ji »  we  obtain 


C?>]  In  =  V„ 


(2-41 ) 


where  I„  and  are  column  vectors.  The  j**-1  element  of  I,-,  is 

b  -9b  b 

Iin-  The  i**'  element  of  Vn  is  given  by 


Vi„  =  <J^E(0,2Mr,)> 


( 2-42 ) 


Furthermore?  CZ**]  is  a  square  matrix  whose  ij^1^  element  is 
given  by 


b  f  b  b 

2‘*  j  il 


?  0)  ds 


( 2-43 ) 


Once  ( Jj  »Mn)  are  chosen?  (2-42)  and  (2-43)  can  be  evaluated. 

-9b 

Thus?  I„  can  be  found  by  solving  (2-41).  Substituting  (2-40) 
into  (2—36)  and  (2-37)?  we  obtain?  in  region  b? 


P*(0?M„)=E(0?2M„)  +  Z  Ij„  E(Jj»0) 

j 


( 2-44 ) 


hf»(0?M„)=H(0,2r%?)  +  Z  IJr,  H(Jj?0) 


(2-45) 


mm 
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5.6  Summary 

In  section  5.4  through  5.5,  we  have  evaluated  all  the 
field  quantities  in  <S-13>— (5— 16) .  Substituting  (5-55)  and 
(5-34)  into  (5-13),  (5-56)  and  (5-35)  into  (5-14),  (5-44) 
into  (5-15),  and  (5-45)  into  (5-16),  we  obtain,  in  region  a. 


Ef"*=  +  E-  =  Elr,*=  +  2  I*  E(  J  j  ,0 ) 

1  •  • 

+  2  V„EE(0,SM„)  +  2  I4„  E  <  J ,  0 )  3 


H*  "*=  +  H-  =  H1  "•=  +  2  I,  H(J.,,0) 

+  2  V„CH(0,3M„)  +  2  Ij„  H(Jj,0)3 

n  J 


and,  in  region  b: 


E*»  -  -  2  V„  CE(0,SM„)  + 


2  E(J*,0)  3 

j 


If  =  -  2  V„  CH(0,5M„>  + 


I  I,n  H ( J j , 0 )  3 
j 


(5-46) 


(5-47) 


(5-48) 


( 5-49 ) 


To  evaluate  V„,  we  must  solve  (5-9)  .  Before  (5-9)  can 
be  solved,  CY*3,  CY*»3  and  Ilnc  have  to  be  evaluated. 
Substituting  (5-35)  into  (5-11)  and  (5-45)  into  (5-15),  we 
obtain 

Ymn  =  Y +  A  V„„  (5-50) 

yI„  =  y2"  A  Y*„  (5-51) 


where 


A  Y, 


» 

Mm  -  H%(  0,SMr,)  ds 
.  A 


-2  I jn  Mm*H*( J^, 

*  J  A 


0)  ds 


( 5-55 ) 


( 5-53 ) 


In  (2— 53),  is  H*(J,,,0)  evaluated  immediately  to 

m 

the  left  of  the  shorted  aperture  on  which  Jj  may  flow.  In 

*►  b  to 

(2-54),  H« ( J j ,0 )  is  H«(Jj,0)  evaluated  immediately  to 

to 

the  right  of  the  shorted  aperture  on  which  J,  may  flow. 

h» 

CYmrt3  is  called  the  half-space  admittance  matrix.  The 
terminology  comes  from  the  fact  that  it  is  the  admittance 
matrix  for  the  case  of  electromagnetic  coupling  through  an 
aperture  in  an  infinitely  large  ground  plane,  which  has  been 
investigated  extensively  in  the  literature,  e.g.,  Bethe  hole 
theory  C81  is  developed  for  the  small  circular  hole  in  the 
ground  plane.  Imagining  the  infinitely  large  ground  plane 
being  shrunk  down  to  a  finite  size  and  bent  over  to  form  a 
cavity,  we  can  view  C  aY*3  and  C  in  (2-50)  and  (2-/1) 

as  modifying  terms.  We  want  to  state  how  Bethe  hole  theory 
should  be  modified  for  an  aperture  in  a  finite  body. 
Substituting  (2-26)  into  (2-10),  we  find  that  the  m**' 
element  of  I1"*1  is  given  by 

I«"C  =  M™,  ds  2  iT*  Mm-HUj,?)  ds 

JA  J  JA  (2-55) 

In  summary,  we  first  shorted  the  aperture  in  the 
conducting  body  with  a  perfectly  conducting  flat  plate  and 
then  put  magnetic  current  sheets  on  both  sides  of  the  plate. 


M  on  one  side  and  -M  on  the  other  to  render  the  tangential 
electric  field  continuous  across  the  aperture.  By  requiring 
the  tangential  magnetic  field  to  be  continuous  across  the 
aperture,  we  found  an  integral  equation  for  M.  Solving  this 
integral  equation  by  the  method  of  moments  C51,  we  obtained 
M  as  the  linear  combination  (2—6)  of  the  magnetic  current 
expansion  functions  Mr,.  In  (2-13>-<2-16) ,  we  expressed  the 
electromagnetic  field  of  the  original  problem  of  Fig.l  as 
the  sum  of  the  field  due  to  the  impressed  sources  and  linear 
combinations  of  the  fields  due  to  the  magnetic  currents  M„, 
all  sources  radiating  in  the  presence  of  the  conducting  body 
with  its  aperture  shorted.  In  turn,  each  of  these  fields  was 
expressed  as  the  sum  of  the  field  due  to  its  source 
radiating  in  homogeneous  space  and  the  field  due  to  the 
electric  current  induced  on  the  body  with  its  aperture 
shorted.  This  electric  current  was  obtained  by  solving  its 
integral  equation  by  the  method  of  moments.  Collecting 
results,  we  were  able  to  express  the  electromagnetic  field 
of  the  original  problem  of  Fig.l  as  the  summations  <2-46)- 
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3 .  ELECTROMAGNETIC  COUPLING  TO  AN  INFINITELY 
LONG  CYLINDER,  TM  CASE 

3 . 1  Remarks  and  Simplifications 

Although  the  formulas  derived  in  Section  £  are  valid 
for  both  two  and  three  dimensional  problems*  it  is  difficult 
to  evaluate  all  the  integrals  in  the  three  dimensional  case. 
Because  of  complication*  we  may  lose  insight  to  the  problem. 
Therefore,  we  shall  make  the  following  simplifications. 

A) .  Everything  is  invariant  in  the  z-direction.  Namely* 
the  conducting  bodies  (three  dimensional  )  become  infinitely 
long  cylinders  and  they  are  completely  specified  by  their  xy 
cross  sections.  The  aperture  becomes  an  infinitely  long 
slot.  The  impressed  sources  ( J*1"*3  ,M*m,:*>  produce  two 
dimensional  fields.  All  field  quantities  in  Section  £  are 
now  only  functions  of  the  coordinates  (x*y).  (The  time 
dependence  exp(jwt)  is  suppressed). 

B)  .  Conductor  Sm  given  in  Fig.  1  is  removed  since  its 
presence  entails  only  minor  modifications. 

C) .  The  conductor  with  an  aperture  which  is  now  a 
slotted  cylinder  has  zero  thickness*  i.e.*  S*.  and  S.  in 
Fig.l  become  indistinguishable  so  that  the  slotted  cylinder 
is  completely  specified  by  one  contour  denoted  by  C  in  the 
xy  plane. 

D)  .  In  Section  3,  we  consider  the  case  where  a  TM  plane 
wave  is  incident  on  the  cylinder.  The  plane  wave  is 
specified  by 
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Elr>e:  =  u,  expC  jk  (  xcosfi^ysinfl1 )  3  (3—1) 

where  u,  is  the  unit  vector  in  the  z  direction*  k=w4n€,  and 
*6*  is  the  polar  angle  of  the  direction  from  which  the  plane 
wave  comes. 

E)  .  As  shown  in  Fig. 5*  the  contour  C  of  the  cylinder  is 
a  chain  of  L-l  straight  line  segments.  The  end  points  of  the 

segment  are  labeled  t  j  and  tj»i  where  t4  is  the  value  of 
t  at  the  beginning  of  the  segment  and  tj*-i  is  the  value  of 
t  at  the  end  of  the  segment.  Measured  from  the  "  beginning  " 
of  C  at  ti=0,  t  is  the  arc  length  along  C. 

F) .  As  shown  in  Fig. 5*  the  contour  ( x=0 *  |  y  |  <.  W  )  of 

the  slot  is  partitioned  into  M  straight  line  segments.  The 
end  points  of  these  segments  are  specified  by  their  values 
of  t.  Here,  t  is  the  arc  length  along  C  extended  onto  the 
contour  of  the  slot.  At  <x=0*  y=-W) *  we  have  t=tL»  and  at 

<x=0,y=W),  we  have  t=ti_*-«  =t,_+EW. 

G) .  The  combination  of  C  and  the  contour  of  the  slot  is 

m  to 

called  C«"= .  Now*  J4  and  J*  are  both  on  C-,=  .  We  choose 


J  j  —  J  j  —  Jj  ( t ) 


<3-5  > 


where 

f u,  t  jit<tj.t 

Ja(t)  ,  j=l *5, . . . ,N  (3-3) 

.  t  0  elsewhere 

where 

N  =  L+M-l  (3-4) 

The  magnetic  current  expansion  function  M„(t)  will  be 


xy  plane  view  of  the  slotted  cylinder.  The  dotted  line 
on  the  y  axis  represents  the  slot.  Csc  is  the  combination 
of  the  contour  C  of  the  cylinder  and  the  contour  (x=0,  |y|  £  W) 
of  the  slot. 
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specified  in  Section  3.3. 

Assuming  that  both  J,»(t)  and  M„(t)  are  tangent  to  Cm<= 
and  independent  of  z,  the  two-dimensional  fields  due  to  the 
electric  current  Jj(t)  on  C**=  and  the  magnetic  current  M„(t) 
on  the  contour  of  the  slot*  both  radiating  in  the  unbounded 
homogeneous  space  with  constitutive  parameters  are 
given  by  C5>Eq .3— £73 » C 1 1 r Eqs . (5 )  and  (6)1 


n 

E( Jj ,0)  =  - 

4 


ti_*(ii 

k  Jj<t’>  Ho  (kj^-^'j)  dt’  + 
0 


tl_»M 


V  dtt-Jjtt’))  «e> 

tJ  ^ — H-‘w 


)  dt' 


0 


(3-5) 


H(  J 


1  r 

j  ,0 )  =  €»  J  j  ( t )  xn  -  -  I 

4j  J 


tu»n 

«*> 

J,(t')  X  vHo  (k|y?-^,J)  dt  ’ 


(3-6) 


Co  Mn  ( t )  xn 


ti_»n 

i  r 

-  Mr,  ( t  ’  )  x  vHo  (  k  I  £—£ *  I  )  dt 

4 j  J  (3-7) 

tL 


1 

H(0,Mr,)= - 

4q 


% 

■I 


t  i_-m 


Mr,(t')  Ho  <  k  |  p-fi  ’  |  )  dt’  + 


ft^.n 

v  d(i*Mr,(t’  )  ) 


k  J  t|_  dt’ 


Ho  <kj^-^’|) 


dt’ 


(3-8) 


where  t  is  the  unit  tangential  vector  on  C**=,  is  the  two- 
dimensional  radius  vector  to  the  point  on  C“e  whose  arc 
length  is  t’,  and  0  is  the  two-dimensional  radius  vector  to 
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the  field  point.  The  field  point  is  the  point  at  which  an 

electromagnetic  field  in  (3—5)— <3— 8)  is  evaluated.  In  (3—5)— 

(3— 8),  k  is  the  wave  number  w -J  y  €  where  w  is  the  angular 

frequency,  rj  is  the  intrinsic  impedance  4h/€  ,  n  is  the  unit 

normal  vector  that  points  outward  from  C®c  as  shown  in  Fig. 5 
«■> 

and  H„  is  the  Hankel  function  of  the  second  kind  of  order 
zero.  If  the  field  point  is  remote  from  C0cr ,  then  €.  is 
zero.  If  the  field  point  clings  to  either  side  of  0®*= ,  then 
€«  is  ±1/2.  Specifically,  €.  is  1/2  if  the  field  point  is  on 
the  region  a  side  of  C*“= ,  and  €.  is  -1/2  if  the  field  point 
is  on  the  region  b  side  of  CSe .  If  C.=±l/2,  then  the  t  that 
appears  in  (3-6)  and  (3-7)  is  the  arc  length  of  the  field 
point  on  C**®. 


3.2  Evaluation  of  Quantities  That  Do  Not  Depend  on 

m 

In  this  section  we  specialize  Z*  *  of  (2-24),  Z,..,  of  (2- 

■  M  •  •=> 

43),  and  V*  of  (2-23)  to  the  case  in  which  both  Jj  and  Jj 
are  the  same  electric  current  J,  of  (3-2). 


a)  Evaluation  of  ZtJ  of  (2-24)  and  Zu  of  (2-43) 
Substituting  (3-2)  into  (2-24)  and  (2-43)  and 
integrating  only  over  the  contour  Cm=  rather  than  a  surface, 
we  obtain 


2,4 


(3-9) 


where 


2,4 


t(.*f1 

Jfc(t)  •  E(J4,0)  dt 

0 


(3-10) 


Substituting  (3-3)  into  (3-5) »  we  obtain 


E(  J j  ,0  )  =  u, 


kq  Pt  <e> 

_  |  H0  <k|/>-£’|>dt’ 

H  J  t  j 


(3-11 ) 


Substitution  of  (3-11)  and  (3-3)  into  (3-10)  yields 


ti»i  t j*i 

kq  f  f 

2*,  = -  dt  dt’  Ho  ( k  J  f>~fL ’  |  > 


(3-13) 


t  *  t  j 

New  variables  of  integration  u  and  u’  are  defined  by 

u  =  3(t-t*-AC*/3)/AC,.  (3-13) 

u’  =  3(t’-tJ-^CJ/S)/^CJ  (3-14) 

where 

j  =  tj*i—  tj  j  j=l*3»....*N  ( 3— 15) 

Substitution  of  (3-13)  and  (3-14)  into  (3-13)  leads  to 


where 

r* t.  =  k AC  i  *  i=l  »3»  ■  .  <  j  N 
r»a(u,u’)  =  kR4J  (r*utt  -  TjU’ijl/S 


(3-16) 


(3-17) 

(3-18) 


In  (3-18) »  Rij  is  the  vector  from  the  midpoint  of  Ca  to  the 
midpoint  of  Ci  where  C».  and  Cj  are  the  straight  line 
segments  of  Cse  whose  end  points  are  (ti»ti4.t)  and 

A  /» 

respectively.  Moreover*  tk  and  are  the  unit 


vectors  tangent  to  Cfc  and  *  respectively.  Here*  t »  points 


from  t»  toward  t4-»,  and  t3  points  from  tj  toward  t  The 

z'  Z' 

vectors  Ri)f  t* *  and  1 3  can  be  expressed  in  terms  of  the 
(x,y)  coordinates  of  the  end  points  of  C*  and  C j . 

If  i  ^  j  *  Zi j  is  found  by  the  two-dimensional  Gaussian 
quadrature  method  C123.  If  i  =  j,  we  recast  (3-16)  as 


n  r*«  1*1  r  1 

!,  %.  =  -  du  du  ’  CH«  (  I  u-u  ’  I  T  ». 

16k  J-l  J-l 

T)  r,-  ri  ri 

j  duj  du’  g(|u-u’jrk/S 


/2)-g( I u  u ’ ir*/2) 3 


(3-19) 


<■> 

where  g(x)  is  the  small  argument  approximation  of  H»  (x) 

C  6  »  p . 462 It  i . e . » 

2 

g ( x )  =  1  -  j - log ( tx/2)  (3-20) 

w 

where  log  denotes  the  natural  logarithm*  logT  is  Euler’s 
constant*  and  t=1.781.  The  first  integral  in  (3-19)  now  has 
no  singularity  at  u=u  ’  and  can  be  evaluated  by  the  two- 
dimensional  Gaussian  quadrature  method.  The  second  integral 
in  (3-19)  is  found  analytically  as  C13*  Eqs.(62)  and  (63)3 


g( ju-u’ |r*/2) 


=  4  +  j4  (  3-21og  <  r\  t/2  )  ) /rr 
Substituting  (3-20)  and  (3-21)  into  (3-19), 


(3-21 ) 

we  obtain 


=  r)  El  +  j(3-21og(r»T/2)  )/ir3/4k 
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ags 


•A< 

*Vri 

•A* 


& 


# 

$ 

i 


v; 


>>x 

sSi 

$ 

A* 

i*?i* 

ft 


&afl 

*3 

*5 

»3 
15 


!*•  *21 
iVY 

Ml 

.♦ii 

•Vrl 

Av 


«V 

•ft 

$ 

i  |V 

•ft: 


r,r,"  fl 


'1  *  < 
du  do  '  CH,_ 

J-l  J-l 


Iju-u’ir^/2)  -  1  -  j— iog(  |u-u’ (T^t/4) 3 

TT 

( 3-se  > 


b)  Evaluation  of  V*  of  (2-23) 

Substituting  <3-l)-(3-3)  into  (2-23) ,  and  integrating 
only  over  the  contour  Cme  rather  than  a  surface,  we  find 
that  the  i**  element  of  V***  is  given  by 


fti^i  jk  (  xcos{6  l+ysin^1 ) 

V*  =  e  dt 

«  1 1 


( 3-23  > 


where  (x,y)  are  the  rectangular  coordinates  of  the  point  t 
on  C-*=.  Performing  the  integration  of  (3-23),  we  obtain 


where 


£c*  sink^i  jkoc* 

V»  - -  e 

k<5* 

«»=[(xt»t+xi  )cos(ll  +  ( y»-^»  +y%  )sin0‘  3/2 
6i  =  C(x1  +  i-Xi  )cos>S1  +  <  y»^»  -y*  )sin^‘  3/2 


( 3-24 ) 


( 3-25 ) 


(3-26) 


Here,  (xt,yt)  are  the  (x,y)  coordinates  of  the  point  t»  on 


Now  that  the  elements  of  Zm  and  V*M  have  been 

— ► 

evaluated,  we  can  solve  (2-22)  for  I“M . 


3 . 3  One  Magnetic  Current  Expansion Function  f g r 

.  the  Narrow  Slot 


MM; 


TO? 


Vil 


In  this  section,  we  assume  that  kw  «  1.  Considering 
that  the  slot  is  narrow,  we  take  only  one  expansion  function 
for  the  magnetic  current  M.  We  choose 


Mi  = 


Uy  -Jl  —  (  y/W) 1 


x=0,  |  y  |  <  W 


elsewhere 


( 3-37 ) 


because  it  is  a  simple  vector  function  that,  for  the  TM 
case,  has  the  correct  direction  and  the  correct  behavior  as 
y  approaches  ±W.  In  the  rest  of  this  section,  formulas 
involving  the  magnetic  current  expansion  functions  are 
specialized  to  the  case  where  M,  of  (3-37)  is  the  only 
magnetic  current  expansion  function. 


a)  Evaluation  of  Vk„  of  (3-33)  and  Vi„  of  (3-43) 
Substituting  (3-3)  into  (3-33)  and  (3-43)  and 
integrating  only  over  Cm=  rather  than  a  surface,  we  find 
that  the  i*1"1  elements  of  Vn  and  V„  are  given  by 


where 


tl_*M 

Vin  =|  Jl (t)  •  E%( 0, 


3fV.)  dt 


( 3-38 ) 


( 3-39 ) 


Referring  to  the  last  paragraph  in  Section  3.5a,  and  using 
(3-7)  with  C6,Eq. (D-15) 3 


<•>  “k(P-/T  >  <e> 

^Ho  <  k  |  |  )  =  - Hx  ( k  |  |  ) 


(3-30) 


substituted  into  it,  we  obtain 


mmm&k 


w 


>sj 

wf 

.•Ml 


■k  r  h  x  (  k  \p-o  > 

—  (Mr,(t’ )x(£-g’ )  )« - - — -  dt’  t  0<  t  <  tL 

!j  J  \r?\ 


E*(0,2M„)= 


(3-31 > 


The  subscript  t  on  both  sides  of  (3—31)  denotes  the 
component  tangent  to  C8c.  Substitution  of"  (3-31)  into  (3- 


29)  yields 


ti_^-ri 


\ 

-  dt Ji ( t ) • 

% 


dt' (M„(t' ) x ( p'  ) 


(kj^'| 


V»„= 


i-1 *2, • • • >L  1 


( 3-32 ) 


i — L 9  L+ 1  * ■ ■ • ■ 


where  N  is  given  by  (3-4). 

Substituting  (3-3)  and  (3-27)  into  (3-32)  and 
introducing  new  variables  of  integration  u  defined  by  (3-13) 
and  u'  defined  by 


u'  =  ( t-tL-W) /W 


(3-33) 


we  obtain 


Vn  = 


jkw^c^ri  ri  Ji-(u' )*  (uB*rl#(u,u'))Hi  ( | r»»(u,u* ) | ) 
-  du  |  du '  ■  ■  - -  -  - - - 

4  J-iJ-i  |rt„(u,u’)| 

9  i=l >2» . . . j L— 1 


( 3-34 ) 


>  i — L  >  L+ 1 > . . . » N 


where 


r».(u,u' )=kR1,+(r»u/2)tl-kWu'uy 


(3-35) 


TO 

? 
f 


•32 

1 


and  R|«  is  the  vector  from  the  origin  at  the  center  of  the 
slot  to  the  midpoint  of  C*. 

Now  that  the  elements  of  Z~ ,  Z** ,  Vx  ,  and  Vi  have 
been  evaluated,  we  can  solve  (2-32)  for  I*  and  <2-41)  for 
Ii .  Because  of  (3-9)  and  (3-28),  the  solution  Ix  to  (2-32) 
is  the  same  as  the  solution  Ix  to  (2-41).  Calling  this 
common  solution  I»,  we  write 

m  b 

I4X  =  Iji  =  Iji  (3-36) 

h> 

b)  Evaluation  of  Ym„  of  (2-52) 

Replacing  the  surface  integral  by  a  line  integral  over 
the  contour  (x=0,  jy|<W)  of  the  shorting  strip,  substituting 
(3-8)  into  (2-52),  and  performing  an  integration  by  parts, 
we  obtain 

ti_>ri  t|_-m 


2knYm„  =  k1 


t|_^-M 


(O 

!  dt  .M„,(t)-  dt  ’  <  t ’ )Ho  <  k  |  p~/3  *  | 

% 


tL 

t 


’  d 

dt  - 

dt 


(  t  -  Mm  ( t )  )  j  dt’ - ( t  •  Mr,  ( t  ’  )  )  Ho  (  k  |  fL—fZ '  J  ) 


( 3-37 ) 


Substituting  (3-27)  into  (3-37),  anticipating  that  t=uv 
when  ti_<t<ti_»M,  and  introducing  new  variables  of  integration 
u  and  u’  defined  by 


u  =  (t-t,_-W)/W 
u’=  (t,-tl_-W)/W 


(3-38) 

(3-39) 


we  obtain 


m. 

ekqYn  =  ( kW )“  Ui  —  Ue 


(3-40) 


where 


fl  _ _  fl  _  ««> 

J du-il-u*  j  du ’ -Jl-(u’  )«  Ho  ( kW  | u-u ’  |  ) 

fl  u  fl  u’ 

du  _  —  du’— -  H„  (kW|u-u’| 

J-l  Jl-u=  J -1  )  = 


(3-41 > 


( 3-42 ) 


We  recast  (3-41)  and  (3-42)  as 


fl  fl  _  <e> 

U»  =  I*  +  duJl-ue  du’ Jl-(u’ )«  Ho  (kWju-u’ | )-g(kW|u-u’ | ) 
J  —  1  4  —  1 


(3-43) 


LU  «  1«  + 


fl  u  fl  u’  ««> 

IcJu — • - Idu’- -  Ho  (kWlu-u’l- 

J-l  4l-u«*  J-l  4l-(u’  )« 


g ( kW  J u-u ’ J ) j 
( 3-44 ) 


where  g  is  given  by  (3-20)  and 


g ( kw | u-u ’ j ) 


=  Jdu-Jl-u®  j'du’-il-(u’  )®  g(kw|u-u’j) 

fl  u  fl  u’ 

du — -  du  ’  — — .  —  g  (  kW  I  u-u  ’  I  ) 

-1  -Jl-u*  J  -1  Jl-7u’  )* 


By  straightforward  integration?  we  have 


(3-45) 


( 3-46 ) 


J-f^ 

|  ubtIi-u' 
fl  u= 

J-l  /l-u« 


du  =  tt/2 


du  =  it/8 


du  =  tt/2 


( 3-47 ) 


(3-48) 


( 3-49 ) 


!5 

B5 

It  is  evident  from  C 14 , Eqs . ( A6a )  and  (A6c)3  that 


f 


1  n 

4 1 — u log  J u-u ’ | du ’  =  -  (u*  -  1/2  -log2)  (3-50) 
-1  2 


It  is  evident  from  C 14, Eq . ( A6b ) 3  that 


fl  u’ 

J-l  -Jl-(u’  )' 


■log  |  u-u’  |  da’  = 


—iru 


(3-51 ) 


Thanks  to  (3—47),  (3—48),  and  (3-50),  I,  of  (3—45)  becomes 

I*  =  ttc/4  +  jw(  l/4-log(TkW/4)  >/2  (3-52) 

Thanks  to  (3-49)  and  (3-51),  I«  of  (3-46)  becomes 

la  —  jTT  (3—53) 

h» 

Now,  2krjYXi  is  given  by  (3-40)  where  U,  and  Ut-  are 
obtained  from  (3-43)  and  (3-44)  in  which  the  explicit 
integrals  are  evaluated  by  the  two-dimensional  Gaussian 
quadrature  method,  and  I»  and  Ie  are  given  by  (3-52)  and  (3- 
53)  . 

a  b 

c)  Evaluation  of  ^Ym„  of  (2-53)  and  ^Ym„  of  (2-54) 

Replacing  the  surface  integrals  in  (2-53)  and  (2-54)  by 


line  integrals  over  the  contour  (x-0,|y|£W)  of  the  shorting 
strip,  we  obtain 
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N  b  to 

=  2  I  jn  Cjm 

j  =  l 


where 


C  jm 


n»<t> 

t._ 

Mm<t) 


I 

V|_ 

f  0 


*  H*(J.,,0)  dt 


H*(Jj,0)  dt 


(3-55) 


(3-56) 


(3-57) 


Substitution  of  (3-9)  into  (3-56)  and  (3-57)  yields 

ti_-*ri 

Cjm  =  J  !3m<t>  ’  dt  (3-58) 

tL 

Cjm  =  -  |  hMt)  •  H%(Jj,0)  dt  (3-59) 

t»_ 

First  substituting  (3-30)  into  (3-6),  then  substituting 
(3—6)  into  (3—58)  and  (3—59),  and  finally  noting  that  J,(t) 
is  on  the  shorting  strip  <ti_  <  t  <  tt_-*.M)  only  when  j  >  L,  we 
obtain 


Cjm  —  Cjm 


tu*«  1 1_ -»■(-( 

r 

|  dt  Mm (  t ) 
t.  0 


jkf  H, 

=  —  dt  Mm ( t )  ■  dt  ’  ( Jj  ( t '  >  x  {prfL7 )  > - - - — 

J 


i  rf  i 

j=l ,2, . . . ,L-1  (3-60) 


Cjm  —  Cjm 


Mm(t) • (Jj(t)xn)  dt 


>  j  — L  j»  L  ^1  <  •  •  »  y  N 
(3-61 > 
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V 

8 


E 


If*  after  interchanging  t  and  t’  in  (3-60),  we  interchange 
the  order  of  integration  in  (3-60),  then  it  is  evident  from 
comparison  with  (3-33)  that 


C  Jm  —  Cjn  —  Vjm  ,j  —  l?2,«..,l—  1 

3 


(3-63) 


Since  Jj(t)  is  given  by  (3-3)  and  Mm  is  Mi  of  (3-37), 
(3-63)  specializes  to 


Cji  =  Cji  —  —  Vji  , j— 1 , 3 , • « « , L— 1 

3 


(3-63) 


where  Vji  is  given  by  (3-34).  Substituting  (3-3)  and  (3-37) 
into  (3-61)  and  introducing  the  new  variable  of  integration 
u  defined  by  (3-3S),  we  obtain 

■J  1-u*  du  ,  j=L ,L+1 , . . . ,N 


•  *.  W  (• 

'4*  *  “  “  J 


( 3-64 ) 


u , 


where 


u4  =  (tj-ti_-W)/W  ,  j=L,L+l  ,  .  .  .  ,N+1 


(3-65) 


Performing  the  integration  in  (3-64),  we  obtain 


a  b  W  —  1  U  j  i 

Cji  =— C j i =  Csin  u  +  u 4 1-u®  1  , j=L ,  L+l , . . . ,N  (3-66) 

4  Uj 


where  the  principal  value  of  sin-1u  is  taken. 

Substitution  of  (3-36),  (3-63),  and  (3-66)  into  (3-54) 
and  (3-55)  gives 


.  1  L— 1  N 

ii  =  —  Z  I ji  Vji  +  Z  Iji  C ji 

a  j=i  j=«- 


(3-67) 
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> 
a 

fcs 


?■ 

s 


»,  1  L-l  N 

4Yix  =  —  Z  Iji  Vji  —  2  Ijx  CJt 

j=L 


( 3-68 ) 


=  -  - 1 


Setting  m=n»l  in  (£-50)  and  (£-51)*  adding  them*  and 
substituting  (3-67)  *nd  (3-60',  we  obtain 


vtx  +  yTx  =  £  YxT  +  2  I,*  VJ4 

j  =  l 


( 3-69 ) 


where  Yu  was  evaluated  in  Section  3.3b*  Ijx  is  the 


element  of  the  column  vector  lx  =  lx  =  I*  that  satisfies 
both  (2-32)  and  (2-41),  and  Vjt  is  given  by  (3-34)  with  i 
replaced  by  j. 


d)  Evaluation  of  Im  of  (2-55) 

Replacing  the  surface  integrals  in  (2-55)  by  line 
integrals  over  the  contour  (  x=0,  jy|  i  W  )  of  the  shorting 
strip*  we  obtain 


t  nc  X  nc  N  bm  m 

1  m  —  C«n  —  2  1  j  C  j  f 

j=l 


(3-70) 


where  Ij  is  the  element  of  the  column  vector  I  that 

m 

satisfies  (2-22),  is  given  by  (3-56),  and 


1 


tl_*M 


X  ne 

Mm  *  H  dt 


(3-71 ) 


tL 


Substitution  of  (3-1)  into  the  Maxwell  equation 

*  nc  j  1 ne 


H 


V  x  E 


( 3-72 ) 


mom 


itnomr?  wn  ruw*- 
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leads  to 

me  1  jk  ( xcos^+ysin**1 ) 

H  =  -  (-U*  sinfi1  -H^y,  cosfi1)  e  (3-73) 

n 

where  u«  is  the  unit  vector  on  x  axis.  Substituting  (3-73) 
into  (3—71)  and  anticipating  that  Mm  will  lie  on  the  y  axis 
and  have  no  x— component *  we  obtain 

t 

cos***j'  jk  ( t-tt.-W)sin|41- 

Cm  = -  ( Mm  -Uy)  e  dt  (3-74) 

n  J 

ti_ 

In  Section  3.3c  *  expressions  (3-61)  and  (3-63)  were  "found 

m 

for  CJn. 

When  h»  of  (3-27)  is  the  only  magnetic  current 

m 

expansion  function*  C,*  is  given  by  (3-63)  and  (3-66)  so 
that  (3-70)  becomes 

Ine  &nc  1  L-l  ■  N  N  ■  x  m 

lx  =  Cx  -  Z  Ij  V,i  -  Z  I,  C,,  (3-75) 

a  j=l  J=L 

where  Cj»  is  given  by  (3-66)  and*  as  obtained  by 
substituting  (3-37)  into  (3-74)  and  introducing  the  new 
variable  of  integration  u  defined  by  (3-38)* 

,„C  w  fl 

Cl  =  -  cos)4l  ■)  1  — u*  cos( kWusin^1 )  du  (3-76) 

q  J-l 

e)  Summary 

When  Mi  of  (3-27)  is  the  only  magnetic  current 

expansion  function,  the  matrix  equation  (3-9)  reduces  to  the 
single  algebraic  equation 


42 


8 


<  Yaa  +  Yaa  >  V»  =  Ij 


( 3-77 ) 


Substituting  <3—69)  and  (3-75)  into  <3-77)  and  solving  for 
Vi,  we  obtain 


Inc  N  v  h 
Ci  —  £  Ij  Cj i 

j=l 


Vi  = 


<3-78) 


r,-  L-l 

2  Ytl  +  2  IJX  VJX 

J  =  1 


-»«M 


where  Ij  is  the  j **  element  of  the  column  vector  I  that 


satisfies  <2-22),  is  given  by  <3-34)  with  i  replaced  by 

•V<i  ^ 

j,  Iji  is  the  j  element  of  the  column  vector  I,  =I»  =  I» 


that  satisfies  both  <2-32)  and  <2-41),  Yaa  is  given  by  <3- 

m  %  nc 


40),  Cji  is  given  by  <3-63)  and  <3-66),  and  Ci  is  given  by 
<  3-76 ) . 

Thus,  when  M*.  <  t )  of  <3-27)  is  the  only  magnetic  current 


expansion  function,  the  constants  Ij,  IJt,  IJt,  and  V*  in 


1  nc  i,  nc 

< 2—46 )-< 2-49 )  can  be  evaluated.  The  fields  E  and  H  are 


given  by  <3-1 )  and  <3-73),  respectively.  The  remaining 

b 

fields  in  <2-46)-<2-49)  are  due  to  the  sources  Jj  and  Jj 
defined  by  <3-2)  and  <3-3)  and  where  is  defined  by 

<3-27).  With  the  meaning  of  2MX  clarified  in  the  last 
paragraph  of  Section  2.5a,  these  fields  are  given  by  <3-5)- 
<3-8).  Having  evaluated  all  quantities  on  the  right-hand 
sides  of  <2-46)-<2-49)  under  the  assumption  that  of  <3- 
27)  is  the  only  magnetic  current  expansion  function,  we  have 
achieved  our  objective  ,  which  was  to  determine  the  total 


fields  (E  +E  *  H  +H  )  and  (E  *  H  )  in  regions  a  and  b. 

3.4  Several  Magnetic  Current  Expansion  Functions  for 

the  Wider  Slot 

In  this  Section*  we  assume  that  the  width  of  the  slot 
is  comparable  to  the  wavelength  so  that  several  magnetic 
current  expansion  functions  are  needed.  The  slot  is 
partitioned  into  at  least  three  straight  line  segments  of 
equal  length  aC.  Otherwise  stated* 


tn-t  -t„  =  aC*  n=L * L+ 1  *  ...  * N  (3-79) 

where  N  is  given  by  (3—4)  in  which  M  >  3.  See  Fig. 5.  The 
magnetic  current  expansion  functions  are  now  defined  by 


r 


4  H(tu.,  -t)/AC)' 

tl_*B  — t 


U> 


M* ( t )  =  { 


*C 


yv 


0 


ti_  i  t  i  ti_*i 

tiL_-.  lit  1  t  L-*-e 

elsewhere 


(3-80) 


Mr*(t  )  = 


t  t|_*n — j 

AC 

t  u.n-*- 1  — t 

AC 

0 


Uy 


tl_+n-  1  — ti  t|_«n 


t  t  —  ti_  +  n.  1 


(3-81  ) 


with  n=a,3, . . . ,M-2 


elsewhere 


m 


jS’il 

$ 


|| 

I? 


I 

s|i 


*vi» 

I 

,*•**»' 


»*i 


ft 

•jSj 

$ 

$ 

I'i'y 
«‘  »• 
s*y 

I'V 

J 

J  M 

m 


«fjJ! 

r 

4 
* 


t— 1|_  t-M—a 


-ti 


MM_1(t)=  -J  1-  ((t-W,)/^)*  u, 


tL<.M 


( 3-82 ) 


elsewhere 


If  M=3,  then  (3-81)  is  to  be  discarded.  The  expansion 
functions  < 3-80 )-( 3-82 )  are  continuous  and  have  the  correct 
behavior  as  t  approaches  t|_  and  as  t  approaches  ti_*-n. 


a)  Evaluation  of  of  (2-33)  and  Vlr,  of  (2-42) 

Equations  ( 3—28 )-( 3—32 )  are  still  valid  because  they 

A  fc» 

were  obtained  without  knowledge  of  Hence,  Vln  and  V*„ 

are  given  by  (3—28)  where  is  given  by  (3-32). 

Substituting  (3-3)  and  ( 3-80 )-( 3-82 )  into  (3-32), 

introducing  the  new  variable  of  integration  defined  by  (3- 
13),  replacing  the  integral  with  respect  to  t’  by  the  sum  of 
the  integrals  over  the  two  straight  line  segments  Ci_  and 

on  which  M„(t’)  exists,  and  introducing  a  new  variable 
of  integration  u’  that  goes  from  -1  to  1  as  t’  goes  from  the 
beginning  to  the  end  of  the  pertinent  straight  line  segment, 
we  obtain,  for  i=l , 2, . . . ,L-1 , 


’  PI  PI 

—  j  k^CjjCi  (Um’C.  h_(u,u’))  <e> 

V»»= -  du  du ' -J4- ( 1-u’  )« - H,(  |rkl_(u,u’  )  |  ) 

16  |rkl_(u,uj)j 

—  1  1 


fl  PI  ( u„  *n  i,t.-*-t(  u ,  U *  )  )  <e  > 

du  du’(l-u’) -  H,  ( 

— i  J— i  | r*  t ( u , u * )  | 


I  r\  *  ( u ,  u  ’ 


f  >  i  > 


( 3-83 ) 


w 


1 


K 


V' 


W 

■%* 


iftr 


J 

i 


vk„= - < 

j  l(b 


pi 
du 
-1 J  -1 


(  y  w  *C  %  t 

du' ( 1+u’ ) - 


-1 (u,u’ )  >  <e> 

Hx< |r1 


jC  1 JL+n — i (u,u’ ) | 


1  (  U  ,  U  '  )  |  } 


1 

du 
-1  J 


PI 

(  UH  *["  x  J  l-»t>  (u,u’ > )  <e> 

du'(l-u') -  Hx  (jr* 


,  ( u  ,  u  ’  )  I  )  3- 


-1 


|C. 


(U.U’ ) 


t  n=2 f 3  , . . . f M— £ 


(3-84) 


j  k^C^Ci 

*1 

Vi  » i"t —  i  —  ” 

du 

16 

\ 

-1  . 

( uN  •  r  x  f 
du' ( 1+u' ) - 

-1  in*. 


— e  (  U  «  U  ’  )  > 


— ( U - U ’ ) 


Hi  <  <u,u')  |  ) 


PI 

du 


1  ( Ux • P x f u»n- i(Uju'))  ce> 

du ’ -J4-(  1+u’ >* - H 1  (  J  r, 

J  —  1  J  —  1  |  Fl  t(Uyj')  | 


,-l(U|U')  j) 

(3-85 


where  Pi4(u»u')  is  given  by  (3-18).  From  (3-38),  we  have 


Vin  =  0 


ji=L,L+l 

^n=l ,2, . . . ,M-1 


( 3-86 ) 


Knowing  the  elements  of  Z* ,  Z** ,  V,-,,  and  V„,  we  can 

4b 

solve  (2—32)  for  Ir,  and  (2-41)  for  I„.  Because  of  (3—9)  and 

-+m 

(3-28),  the  solution  I„  to  (2-32)  is  the  same  as  the 
solution  I,-,  to  (2—41).  Calling  this  common  solution  I„,  we 
wr  i  te 


I  in  —  I  Jr>  —  I 


b>  Evaluation  of  Ymn  of  (2-52) 


(3-87) 


Equation  (3—37)  is  still  valid  because  it  was 
obtained  without  knowledge  of  ft*  and  Substituting  (3- 


80)-(3— SE)  into  (3—37),  replacing  the  integral  with  respect 
to  t'  by  the  sum  of  the  integrals  over  the  two  straight  line 
segments  and  on  which  M„(t’>  exists?  .introducing 
a  new  variable  of  integration  v’  that  goes  from  0  to  1  as  t’ 
goes  from  one  end  to  the  other  end  of  the  straight  line 
segment?  and  dealing  similarly  with  the  integral  with 
respect  to  t  to  introduce  a  corresponding  variable  of 
integration  v?  we  find  that  the  first  row  of  Y*'”*  is  given  by 


2knY**  =(k^C)=  dvJl-v«  dv’Jl-v'=Ho  (k<iC|v-v'|) 

J  0  J  0  '  ' 


fl  fl  _  <«> 

+2  dv  v  dv’J  l-v’~“  H0  (kAC|  v-v’-l  |  ) 

,  0  J  0 

fl  fl  <«> 

+  dv  v  dv’v’  H0  (kAC|v-v’|> 

J  0  J  0 

*1  v  fl  v’  <ls> 

-  dv -  dv' -  H0  (  k |  v— v ’  |  ) 

J  0  Jl-v=  J  0  Jl-v’  = 

‘1  fl  v'  « e  >  fl  fl 

2  dv  dv' - -  -H0  < k£C I v-v ' -1 |  )  -  dv  dv’  Ho  (k£ 

0  J0  4l-v’«  J0  J  0 


( 3-88 ) 


Cl v-v’ I ) 


2kr)Ylr,=  (kAC> 


fl  _ fl 

dv4l-va  dv' 

J  0  J  0 


v '  C Ho  <  kMZ  )  v+v  '  +n-2  |  )  +H0  (  k |  v-v  ’  +n  j) 


fl  fl  <s>  <*=> 

+  dv  v  dv’v’CHo  (k^C| v+v'+n-3| )+Ho  ( k^C | v-v ’ +n-l j >  ] 

J  0  J  0 

fl  v  fl 

- | dv- - j  dv ' CHp  ( k^C | v+v  7  +n-2  j  ) -H0  ( kAC I v-v ' +n I >  1 


■<ry 


ml 


Since  t  and  t ’  can  be  interchanged  without  affecting 


\p-fi.*  |  in  (3-37),  is  symmetric. 

ri«  ha 

=  Y„m  (3-92) 

h« 

If  both  Mm  and  are  reflected  about  the  xz  plane,  then  Ym„ 
does  not  change.  Since  the  reflection  of  Mm  is  Mm_„,  and  that 
of  M„  is  Mm— „,  it  is  evident  that 


ha 

Ym — m  »M — n 


( 3-93 ) 


From  (3—92)  and  <3— 93),  we  obtain 


Y„x  =  Yjr, 


Yn-l,n  —  Ytfn-n 


n=l , 2 , . . . , M— 1 


(3-94) 


Yr.  ,  r-l—  l  —  Yi(M_n  J 


Therefore,  the  first  column,  the  last  row,  and  the  last 
column  of  Yh*  can  all  be  generated  from  the  first  row  of  V*" 
given  by  ( 3—88 )-( 3-90) .  It  is  not  necessary  to  compute  all 
the  internal  elements  of  Y*’*  because,  as  given  by  (3-91), 
they  depend  only  on  ] m-n j  rather  than  on  m  and  n 
individual ly . 

Each  integral  in  ( 3-88 )-( 3-91 )  whose  integrand  does  not 
have  any  singularity  interior  to  the  region  of  integration 
is  evaluated  by  the  method  of  two-dimensional  Gaussian 
quadrature.  Each  integral  that  has  a  singularity  interior  to 
the  region  of  integration  is  treated  in  Appendix 


A  by 


"subtracting  out"  the  singularity  before  applying  the  method 
of  two-dimensional  Gaussian  quadrature. 


c)  Evaluation  of  of  (2-53)  and  ^Ym„  of  (2-54) 

What  was  done  in  the  first  two  paragraphs  of  Section 
3.3c  is  still  valid  because  it  was  accomplished  without 
knowledge  of  the  magnetic  current  expansion  functions.  In 

4i  b 

those  two  paragraphs,  we  found  that  ^Ymr,  and  ^Ym„  are  given 

A  fc» 

by  (3—54)  and  (3—55)  where  Cjm  and  CJm  are  given  by  (3-62) 
and  (3-61). 

Since  Jj(t)  is  given  by  (3-3)  and  M„,  is  given  by  ( 3- 
80 ) - ( 3-82 ) ,  ( 3-62 )  spec i a 1 i zes  to 


«  1  !j— 1,2,...,L  1 

(-jm  —  Cjm  =  Y  j  m  j/ 

a  m=l ,2, . . . ,M— 1 


(3-95) 


where  is  given  by  ( 3-83 )—( 3-85 ) .  Substituting  (3-3)  and 
( 3—80 )—( 3-82 )  into  (3-61)  and  evaluating  the  resulting 
integrals,  we  obtain 


=  -Ct_i  = 


Si_ i.i  — 


— 1  •  m  Ct_xn — 1  •  m  “ 


a  an  - 


( 3—96 ) 


(3-97) 


m— 2 , 3 ■ 


run* > 


M 


m 


Bjn— 1  “ 

4 


Ci_-*-ri—  x  f  ri—  x  —  i  (M-i  — 


( 3-98) 


The  rest  of  the  C~’s  and  Cto’s*  those  that  do  not  appear  in 
( 3-95 )-( 3-98) »  are  zero. 

Substitution  of  (3-87),  and  <3-95>-<3-98)  into  (3-54) 
and  ( 3-55 )  g i ves 


where 


a  1  L— 1 

AYmn  —  —  21  1  jn  Vjm  +  /\Y  r 

2  j=l 

„  1  L— 1 

^Ymn  —  21  I  jn  V  r 

2  j-1 


^  AC  tr 

*— (  —  it _ n  ^  r  n  ) 

4  2 


Yntn  —  A  "  (  2  L.  -*-m —  1  |  n  ^  1  r  ri) 

4 


(  I u-*-n — e r  n  +  — —  rn) 

4  2 


(3-99) 


(3-100) 


I  ffl - 2  f  3  f  m  m  m  ,  fl  2 


(3-101 ) 


, m=M-l 


Substituting  (3-99)  and  (3-100)  into  the  sum  of  (2-50)  and 
(2— 51 ) »  we  obtain 


» mn  T  » mn 


L-l 

2  Ym„  +  Z  I  jr,  V„ 


(3-102) 


where  Y«n  was  evaluated  in  Section  3.4b,  I  Jn  is  the  j'”-' 
element  of  the  column  vector  I„  =  I„  that  satisfies  both  (2- 


SOfivW 


51 


K 


f 


I 

6 


*> 


32)  and  (2-41),  and  VJm  is  given  by  (3-83)-(3-85)  with  i 
replaced  by  j. 


i  nc 

d)  Evaluation  of  Im  of  (2-55) 


What  was  done  in  the  first  paragraph  of  Section  3.3d 
is  still  valid  because  it  was  accomplished  without  knowledge 
of  the  magnetic  current  expansion  functions.  In  that 

i  nc  •  h 

paragraph,  we  found  that  Im  is  given  by  (3-70)  in  which  I  4 
is  the  element  of  the  column  vector  I***  that  satisfies 

m  1  nc 

(2-22),  CJm  is  given  by  (3-56),  and  Cm  is  given  by  (3-74) 


When  the  magnetic  current  expansion  functions  are  given 
by  ( 3-80 )-( 3-82 ) ,  the  non-zero  C*«.’s  are  given  by  (3-95)-(3- 
98)  so  that  (3-70)  becomes 


„  AC 

TT  •  x  sex 

r _ 

(  -  Ii_  +  Ix_~x  ) 

,m=l 

4 

2 

(3-103) 

l.  nc  1  nc  1 

L— 1 

AC 

Im  s  Cm  “  “ 

Z  Ij  V  jm— < 

\  ~ 

(  I  l_-s-m —  1  ^  I  l_  -*-m  ) 

,  m=2 , 3 ,  . .  .  ,  M— 2 

5 

j  =  l 

4 

1 

!  AC 

x 

i 

i 

i 

|  - 

(  Ii_*n-a  +  ”  1 i_ 

X  )  ,m=M-l 

:  4 

2 

Substituting  ( 3— 80)— (3— 82)  into  (3—74),  we  obtain 


Cr  = 


ACcosj*1 


n 


ri  _ 

i  J  l-v® 
,  0 


-j#lV 

e  dv 


ri  -j**v 

v  e  dv  3 

J0 


(3-104) 


•■.‘XI 


si 


^CCOS#1  j  ( m#  i+Pe )  |  -j#»fl  j#iV 


j#i  fl  ~j#iV 


(3-105) 


m»2 1 3  $ ■ • • t M— 5 


»„«=  ^CCOS#1  j( (M~l )#!+#=)  fl _  j#iV 

Cn-1  *  -  e  C  j  l-va  e  dv 

rj  J  0 


_j#X  fl  j  01V 


dv  3 


(3-106) 


where 


# 1  *  k^C  sin#1 


#«  -  -  KW  sin#1 


(3-107) 


(3-108) 


If  sin#1  #=  0,  then  #*#=  0  and*  thanks  to  E15, Formula  567.1.3, 
( 3— 1 04 )  —  <  3— 1 06 )  become 


me  ^Ccos#1  j(#i+#e)f  PI  _  —  j  #  i  v 

Ci  =  -  e  -i  l-ve  e  dv 

J  0 


1-cos#!  j(#i-sin#i) 


(3-109) 


me  2^C  cos#1  j(m#i+#e) 


Cm  = 


r]  #i' 


C  1— cos#i  3  , m— 2 , 3 , ... »M— 3 


(3-110) 


Cr,-i  - 


^Ccos# 1  j ( (M— 1 )# 


i+#«) r  ri 

I  ■> 

j0 


j#»v 


l-cos#i  -  j(#i  -sin#») 


(3-111 ) 


VLVWJV 


If  sin^*  =  0,  then  =#e  =0  and  (3-104)-(3-106)  reduce  to 


ft*.  = 


AC  cos^1 

-  < 

2q 

AC  cos** 


n 


TT 

1  +  -  > 

2 


,  ws2 1 3 , 


t  M-H 


(3-112) 


(3-113) 


1  nc  1  nc 

Cn- i  =  Ci 


(3-114) 


l  nc  •  k 

Now,  Im  of  (2-55)  is  given  by  (3-103)  in  which  I ^  is 


the  j**  element  of  the  column  vector  I-**  that  satisfies  (2- 

1  nc 

22),  V.,m  is  given  by  ( 3-83 )-( 3-85 )  ,  and  Cm  is  given  by 
either  ( 3-109)-(3-l 1 1 )  or  (3-112)-(3-114) . 


e )  Summary 

Substituting  (3-102)  and  (3-103)  into  (2-9),  we  can 
solve  (2-9)  for  V,  whose  nt,r'  element  V„  appears  in  (2-46)- 

•  x 

(2-49).  Of  the  remaining  constants  in  ( 2-46 )-( 2-49 ) ,  Ij  is 
the  element  of  the  column  vector  I  that  satisfies  (2- 

—  b 

22),  and  Ij„  and  are  given  by  (3-87)  in  which  IJr.  is  the 

->to 

element  of  the  column  vector  I  =1=1  that  satisfies 
both  (2-32)  and  (2-41). 

Thus,  when  the  magnetic  current  expansion  functions  are 

•  m  «  b 

given  by  < 3-80) -( 3-82) ,  the  constants  V„,  Ij,  and 

in  ( 2-46 )-( 2-49 )  can  be  evaluated.  The  fields  ElrM=  and 
are  given  by  (3-1)  and  (3-73),  respectively.  The  remaining 

m  b 

fields  in  ( 2-46 )-( 2-49 )  are  due  to  the  sources  Jj  and  J, 
defined  by  (3-2)  and  (3-3)  and  2M„  where  is  defined  by 
< 3-80 )-( 3-82 ) .  With  the  meaning  of  2M,-,  clarified  in  the  last 


4 .  NUMERICAL  RESULTS  AND  DISCUSSION 


Programs  have  been  written  in  FORTRAN  for  both  a  narrow 
slot  and  a  wide  slot.  More  general  than  the  theory  presented 
in  Sections  3.3  and  3.4,  these  programs  apply  when  the  wall 
of  the  cavity  is  finitely  thick  and  when  an  additional 
conductor  is  either  inside  or  outside  the  cavity.  A  report 
on  the  usage  of  the  programs  is  in  preparation.  Since  the 
magnetic  current  (or  the  aperture  field)  is  our  main 
interest,  it  has  been  computed  for  different  widths  of  the 
aperture,  different  sizes  of  the  cavity  and  different  angles 
of  incidence.  Some  far  field  patterns  are  also  computed. 

4. 1  Remarks  and  Definitions 

The  numerical  results  presented  in  this  report  are  all 
for  the  case  where  the  uniform  TM  plane  wave  given  in  (3-1) 
impinges  on  the  aperture  with  the  incident  angle  being 
defined  in  Fig. 5.  We  decided  to  plot  the  magnitude  of  the 
magnetic  current  versus  the  position  in  the  aperture  since 
our  computations  show  that  the  phase  of  the  magnetic  current 
changes  little  in  the  aperture.  All  numerical  integrations 
are  performed  by  10  point  Gaussian  quadrature. 

The  name  pseudo-image  method _ 1 _ is  attached  to  the 

solution  developed  in  Section  3.3,  where  only  one  special 
expansion  function  is  used  for  the  magnetic  current.  The 
name  pseudo— image  method  S _ is  attach ed  to  the  solution 


developed  in  Section  3.4, 


where  no  less  than  two  expansion 


functions  are  chosen  for  the  magnetic  current.  The  Fourier 
Series  method  is  the  method  presented  in  Cl 3. Appendix  B3  for 
a  circular  cylindrical  cavity  with  very  narrow  slot.  The 

scattering  method _ is  the  method  used  in  E163  and  summarized 

in  Appendix  B,  where  only  an  electric  current  is  solved  for. 
In  the  scattering  method,  a  matrix  equation  is  extracted 
from  the  electric  field  integral  equation  either  by 
Galerkin's  method  C5, Section  1—33,  or  by  point  matching.  In 
Galerkin’s  method,  the  symmetric  product  of  the  integral 
equation  is  taken  with  each  expansion  function.  In  point 
matching  C5, Section  1-43,  the  integral  equation  is  enforced 

at  discrete  points.  The  non-pseudo- image  method _ is  the 

method  where  the  pseudo-image  introduced  in  Section  2.5  is 
eliminated.  This  method  is  summarized  in  C13,  Appendix  A3. 
Finally,  if  the  plots  are  for  the  pseudo-image  method,  then 
the  magnetic  current  expansion  functions  used  are  triangles 
only  (i.e.,  M»(t)  of  <3-80)  and  t ( t )  of  (3-82)  are 
replaced  by  triangles)  unless  we  state  otherwise.  The  reason 
for  this  is  explained  in  Section  4.4. 

4.2  Validity  of  Results 

We  have  to  check  whether  we  obtain  correct  and 
accurate  results.  This  is  not  only  an  important  task  but 
also  a  difficult  one  since  no  exact  solution  is  available  to 
compare  with.  That  is  why  we  started  with  the  simplest  case: 
a  circular  cylindrical  shell  of  zero  thickness.  This  shell 
was  treated  by  the  Fourier  series  method  in  C133  where  the 
magnetic  current  distribution  (3-27)  was  on  the  arc  in 


Fig. 6  instead  of  on  the  plane  strip.  For  ease  of  comparison 
and  simplicity*  most  of  our  computations  are  made  for 
circular  cylinders  and  zero  thickness  of  the  cavity  wall 
even  though  the  programs  work  for  thick  cylinders  of 
arbitrary  cross  section.  The  comparison  is  made  as  follows: 

It  is  very  interesting  to  note  the  following  from 
Figs. 7— 16: 

a)  The  results  obtained  by  the  different  methods  are  of 
the  same  order  of  magnitude  (Fig. 7).  The  pseudo— image 
methods  and  the  Fourier  series  method  give  especially  close 
results.  For  instance*  at  the  center  of  the  aperture  the 
magnetic  current  is  obtained  as  0.15234  /-103. 13°  by  the 
pseudo-image  method  1*  0.15384  /-103.490  by  the  pseudo-image 
method  8*  and  0.1535  /-105.970  by  the  Fourier  series  method. 

b)  The  scattering  method  yields  larger  results  than 
other  methods.  Scattering  (Galerkin’s  method  )  yields  the 
second  largest  and  Scattering  (point  matching)  yields  the 
largest.  (see  Figs. 3-10) .  The  results  obtained  by  the 
scattering  methods  monotonical ly  approach  the  results 
obtained  by  the  pseudo-image  method  as  N  increases.  This  is 
true  for  both  a  large  aperture  <e.g.»  ^*,=30°)  and  a  small 
aperture  (e.g..  ^o=50 ) .  (see  Figs. 11-16).  Furthermore. 
Scattering  (Galerkin’s  method)  approaches  faster  than 
Scattering  (point  matching).  (see  Fig. 13).  The  magnetic 
currents  of  both  scattering  methods  are  usual  lv  within  l&'A 
of  each  other.  Scattering  (Galerkin’s  method;  is  more 
accurate  since  Galerkin’s  method  tests  the  equation  over  an 


An  infinitely  long  circular  cylinder  with  a  slot.  k=  (Ov'UC  is  the 

wavenumber  and  a  is  the  radius  of  the  circular  cylinder. 

-f  +  $  <d)<Tr-(t>  defines  the  conducting  oart  of  the 

o  —  —  o 

cylinder.  The  edges  of  the  shell  are  at  (x  =  0,  y  =  ±  a  sin  ) 
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interval  instead  of  at  one  point  as  point  matching  does 
C5, Sections  1-3  and  1-43. 

As  a  further  check.  we  compare  our  results  with  those 
obtained  in  C163  and  C173  for  two  special  cases.  At  the 
center  of  the  aperture,  the  magnitudes  of  the  electric  field 
are  obtained  as:  case  i)  0.31  by  C163  and  C173,  0.335  by 
Scattering  (point  matching)  with  N=31 .  0.33  by  Scattering 
(Galerkin’s  method)  with  N=31  ,  and  0.19  by  the  pseudo-image 
method  with  M=15  and  L=18>  where  $4l=0°,  560=10°, and  ka=l  . 
Case  ii)  0.75  by  C163  and  C173,  0.54  by  Scattering 
(Galer kin’s  method)  with  N=31 ,  and  0.53  by  the  pseudo- image 
method  with  M=15  and  L=18»  where  01=0°,  j6o=30°,and  ka=l  . 
f*1  ,  *o»  and  ka  are  defined  in  Fig. 6. 

In  the  previous  paragraph,  the  numbers  obtained  in  C163 
and  E173  are  the  largest.  The  reasons  are  twofold.  First  of 
all,  C163  and  C173  evaluate  the  electric  field  on  the  arc 


I 


:! 

: 


part  of  the  cavity  (Fig. 6  ),  whereas  we  evaluate  the 
electric  field  on  the  plane  strip  connecting  the  cavity 
edges.  Ours  should  be  smaller  since  we  are  further  inside 
the  cavity.  However,  This  should  not  yield  a  big  difference 
when  the  aperture  is  small.  Secondly,  in  E163  and  C173  point 
matching  is  used  to  apply  the  moment  method,  and  the 
integrals  are  approximated  by  the  interval  of  integration 
times  the  integrand  sampled  at  the  center  of  the  interval 


with  special  treatment  of  the  singularity  of  the  integrand. 
Gaussian  quadrature  yields  results  that  are  more  accurate 
than  those  of  any  other  technique  for  well  behaved  functions 


Ciai-  Hence,  our  results  should  be  more  accurate  than 
theirs . 

On  the  whole, the  comparisons  and  discussion  presented 
above  lead  to  the  conclusion  that  the  results  we  obtained 
are  correct  and  that  the  pseudo-image  method  yields  very 
accurate  results. 

4 . 3  Usefulness  of  the  Pseudo-image 

On  the  surface,  the  introduction  of  the  pseudo-image  in 
Section  2.5  seems  a  little  artificial.  However,  as  shown  in 
Figs. 17— 19,  the  pseudo-image  method  does  give  better  results 
than  the  non-pseudo-image  method.  First  of  all, the  non- 
pseudo-image  method  yields  unexpected  overshoots  near  the 
edges.  We  say  that  they  are  unexpected  overshoots  since  no 
other  method  mentioned  in  Section  4.1  predicts  them. 
Secondly,  the  non— pseudo-image  method  gives  larger  amplitude 
(within  10%)  but  essentially  the  same  phase  for  the  magnetic 
current  (within  1 %  ). 

From  the  results  discussed  in  Section  4.2  and  those 
shown  in  Figs. 17-19,  we  conclude  that  the  pseudo— image 
method  gives  a  more  accurate  magnetic  current  than  the  non- 
pseudo-image  method. 

4 . 4  Edge  Conditions 

In  Section  3.4,  we  have  chosen  the  two  special 
expansion  functions  Mt(t)  and  l,lM_1(t>  to  obtain  the  proper 
edge  behavior  C141  because  we  initiallv  believed  that  this 
would  yield  more  accurate  results.  Alternatively,  we  could 


use  only  triangle  functions.  For  comparison  we  carry  out 
both  solutions.  Figs.  20-22  show  that  the  results  of  these 
two  methods  are  very  close  to  each  other?  even  near  the 
edges.  Not  much  advantage  is  gained  by  satisfying  the  edge 
conditions  in  the  examples  we  computed.  The  addition  of  the 
special  expansion  functions  results  in  much  more  complicated 
equations?  and  more  effort  and  computer  time  are  needed. 

4.5  Speed  of  Convergence 

Given  a  certain  size  of  aperture?  how  many  expansion 
functions  should  be  used?  In  other  words?  how  fast  do  the 
results  converge  ?  It  can  be  seen  from  Figs. 23— 26  that  M=10 
suffices  to  obtain  convergence  for  a  small  aperture  and  M=16 
suffices  for  a  large  aperture.  The  larger  M  is?  the  smoother 
the  plots  are.  Hence?  as  a  rule  of  thumb?  M  should  be  large 
enough  to  obtain  a  smooth  curve. 

4.6  Other  Numerical  Examples 

Fig. 27  shows  that  the  magnitude  of  the  aperture  field 
increases  with  the  width  of  the  aperture.  Figures  28  and  29 
show  that  oblique  incidence  causes  noticeable  asymmetry  of 
the  aperture  field  only  when  the  aperture  is  large.  Figs. 30 
and  31  give  results  for  the  case  where  the  cavity  wall  is 
finitely  thick  (see  Fig. 32).  It  is  interesting  to  note  that 
the  aperture  field  increases  with  the  wall  thickness. 

The  far  scattered  field  is  obtained  by  replacing  the 
Hankel  function  by  its  large  argument  approximation  C6? 
eq.(D-13)l.  Although  the  aperture  fields  obtained  by  the 


different  methods  mentioned  in  Section  4.1  are  d i f f erent , the 
far  fields  are  nearly  the  same.  Their  difference  is  less 
than  1 V,.  This  is  expected  since  the  aperture  field  only 
affects  the  electric  current  in  the  vicinity  of  the 
aperture.  However,  the  field  inside  the  cavity  is 
proportional  to  the  aperture  field.  Examples  of  far  field 
scattering  patterns  are  given  in  Fig. 33  and  Fig. 34.  If  the 
angle  of  incidence  f«$*  increases  by  •  then  the  far  field 
scattering  pattern  rotates  through  the  angle  <5*4*  as  does  the 
scattering  pattern  of  an  infinitely  long  complete  circular 
cylinder.  Again,  this  shows  that  the  small  aperture  has 
little  effect  on  the  far  field  scattering  pattern. 

4.7  Concluding  Remarks 

In  this  report,  a  new  method  called  the  pseudo-image 
method,  incorporated  in  the  generalized  network  formulation 
for  aperture  problems,  is  developed  to  accurately  determine 
the  field  inside  and  outside  a  conducting  cavity  with  a 
small  aperture.  The  theory  can  be  applied  to  both  two  and 
three  dimensional  cavities  of  arbitrary  shape. 

A  number  of  computations  have  been  made  for  an 
infinitely  long  slot  in  a  perfectly  conducting  cylindrical 
surface  illuminated  by  a  uniform  TM  plane  wave.  The  programs 
developed  apply  to  a  cavity  with  finite  thickness  and  an 
additional  conductor  may  be  either  inside  or  outside  the 
cavity.  With  some  modifications,  the  main  program  can  apply 
to  a  multi-conductor  system  as  well.  As  further  work,  we 
could  solve  the  problem  of  an  infinitely  long  cylinder  with 


a  slot  illuminated  by  a  uniform  TE  plane  wave  instead  of  a 
uniform  TM  plane  wave. 

Final ly,  we  should  point  out  that,  although  the  pseudo¬ 
image  method  works  well  for  a  non— resonant  cavity,  it  fails 
when  the  cavity  is  resonant.  There  two  reasons.  One  is  that 

Sc 

the  resonant  electric  current  on  the  conducting  surface  S„ 

produces  tangential  electric  field  on  this  surface  so  that 
a 

‘Urt  can  not  be  uniquely  determined  by  using  the  electric 
field  integral  equation  only  C18,  Section  S3.  The  other  is 
that  the  field  inside  the  cavity  would  go  to  infinity  if  the 
magnetic  current  on  the  right  hand  side  of  the  shorted 

aperture  in  Fig. 4  excited  the  resonant  mode  of  the  cavity.  A 
special  technique  has  to  be  developed  to  treat  the  resonant 
cavity . 
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0  Pseudo-image  method  (M=16,  L=33) 

X  Scattering  (Galerkin's  method)  N=31 
f Scattering  (point  matching)  N=31 
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The  magnitude  of  the  magnetic  current  in  the 
aperture  for  =  180°,  =  10°,  ka  =  1 . 
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Fig.  12.  The  magnitude  of  the  magnetic  current  in  the  aperture  for 


l)1  -  180°,  <*>  =  30°,  ka  «  1. 
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Fig.  18.  The  magnitude  of  the  magnetic  current  in  the  aperture 


for  =  180°,  i>  =  i0° *  ka  =  1 . 
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Fig.  20.  The  magnitude  of  the  magnetic  current  in  the  aperture  obtained  by  the 

pseudo  image  method  with  and  without  the  satisfaction  of  the  edge  condition 

i1  =  180°,  i  =  1.25°,ka  =  2. 
o 


0  Triangles  only 
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Fig.  22  . 
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The  magnitude  of  the  magnetic  current  in  the  aperture  obtained  by 
the  pseudo-image  method  with  and  without  the  edge  condition 

included  for  <i>^  =  180°,  t>  ~  30°,  ka  =  1 . 
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Fig.  25.  The  Speed  of  convergence  of  the  magnitude  of  the  magnetic  current 
in  the  aperture  for  the  pseudo-image  method  with  L  *  18,  4)1  *  180 
<t>  *  15°,  ka  ■  1. 
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Fig.  26.  The  Speed  of  convergence  of  the  magnitude  .  F  the  magnetic  current 

in  the  aperture  for  the  pseudo-image  method  iwth  L  *  18,  4)1  *  180°, 
$  =  15°,  ka  =  3. 


The  magnitude  of  the  magnetic  current  in  the  aperture  obtained 

by  the  pseudo-image  method  for  different  incident  angles,  where 

1)  a  5°,  ka  3  1. 
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Fig.  29.  The  magnitude  of  the  magnetic  current  in  the  aperture  obtained  by 

the  pseudo-image  method  for  the  different  incident  angles,  where 

1)  *  30°,  ka  =  1. 

o 


Fig.  30. 


The  magnitude  of  the  magnetic  current  in  the  aperture  obtained  by  the 
pseudo-image  method  for  a  cylinder  with  finite  thickness,  <f>-*-  =  180°, 


=  5' 


ka:  =  1. 
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Fig.  31.  The  magnitude  of  the  magnetic  current  in  the  aperture 

obtained  by  the  pseudo-image  method  for  a  cylinder  with 

finite  thickness,  =  180°,  1)  =1 5°,  kan  =  1. 
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Fig.  33.  Far  field  scattering  pattern  for 
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APPENDIX  A 


In  this  Appendix^  we  express  2krj  Ymn  of  (3-88)-(3- 

91)  in  terms  of  integrals  whose  integrands  have  no 

singularity  interior  to  the  region  of  integration.  Whenever 

<  • » 

the  argument  of  H0  is  zero  at  a  point  interior  to  the 
region  of  integration  we  subtract  out  the  small  argument 

<  m  > 

approximation  of  Ho  and  integrate  this  approximation 
analytical ly . 

Of  all  the  arguments  of  the  Hankel  functions  in  (3— 88), 
only  kAC|v-v*j  can  be  zero  at  a  point  interior  to  the  region 
of  integration.  Accordingly,  we  recast  (3-88)  as 


2krjY»  »  =  <  k^C  )> 


[  fdv  [J 

J0  J0 

- 


civ*  ( 4 1— v*  4l-v y  m  +VV *  ) 
0 


(Ho  ( k^C | v-v * J )  -  g ( k^C | v-v’ | )  > 

fl  fl  _ _  <«> 

!  dv  v  dv’4l-v’e  Ho  ( k^C) v-v’-l | )  +1, 

J0  J  0 


( A— 1 ) 


‘1  PI  v  v' 

—  dv  dv’ ( 

0  0  Jl-V*  -/ 1  —v  ’ 1 


+1)(H0  (k^Cjv-v’j)  -  g(k^C| v-v’ | ) ) 


PI  PI  v’ 

l  dv  dv’ - Ho  (k^CJ  v-v’-l  j  )  - 

J  0  J  0  Jl-v’  = 


where  g  is  given  by  (3-20)  and 


PI  PI 

. »  =  dv  dv 
J0  J0 


( '1 1  — v*  -J 1  -v  ’  *  +■  vv  ’  )  g  (  k^C  j  v-v  ’  |  ) 


( A— 2 ) 
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-  [iv  fi 

J  0  J0 


V  V  7 

dv7  ( - 

0  Jl-vB  -il-V  : 


+  l)g(kiic|v-v7  |  ) 


( A— 3 ) 


The  indentity 


fdv  fi 

j0  J  0 


dv 7  -/i -v  ’ 
0 


log  J  v-v 7  J  = 


( A— 4 ) 


PI 


Er 


dv-Jl— v*  j  dv  * 4 1-v 7 *log  J  v-v 7 


-I 


PI 


dvT 1 -vfc 
0 


dvJ  -J1— v  7  *1  og  I  v+v 7 
0 


can  be  verified  by  expressing  each  integral  over  (-1,1)  on 
the  right-hand  side  of  (A-4)  as  the  sum  of  the  integral  over 
(-1,0)  and  the  integral  over  (0,1)  and  changing  the  variable 
of  integration  so  as  to  replace  the  integral  over  (-1,0)  by 
an  integral  over  (0,1).  Using  (3-47),  (3-48),  and  (3-50),  we 
obtain 


[: 


PI 


dvil-v1 

-1 


IT * 


d  v  *  -J  1  -v  7  *  log  I  v-v' I  - - (1+  41og2  ) 

J-l  16 


( A— 5) 


Substitution  of  (A— 5)  into  (A-4)  gives 


J 


1  _  PI 

dv-Jl -vB  dv  7  -Jl-v 7  *  log  I  v-v : 
0  0 


PI  . _  PI  _ _ _ 

-  -  (  1  +  41og2  )  -  dv  -JT-V*’  dv 7 -Jl-v 5 *  log  I  v+v 7 

32  J  0  J  0 


( A— 6 ) 


Using  C15, formulas  610.,  610.1.,  610.2.,  and  610.3.3 

and  working  diligently,  we  obtain 

7 


‘1  PI 

dv  v  d 
0  0 


dv7  v7  log  I v-v 7  I  - - 

0  16 


( A— 7 ) 


It  is  easy  to  show  that 


Substituting  (3-20)  into  <A-2>  and  using  <A-6)-(A-B),  we 
obtain 


ti*  1  it  1  n  it  7 

I x  - -  +  -  -  j  (  (-  +  - )  1  ogb  -  — -  log2 - - 7-  ) 

16  4  8  2n  4  16  8w 

j2  fl  _  PI  _ 

+  "  dv-/ 1  —v*  1  d v  7  -fl- v 7 “  1  og  I  v+v 7  [  ( A— 9  > 

TT  J  0  J  0 


where  t 

b  =  -  k^C  (A- 10) 

2 

Seeking  to  evaluate  X*  of  (A-3),  we  consider 


PI  v  PI  v' 

U*  =  }  dv  -  ■  I  dv7 -  log|v-v7|  (A-ll) 

J0  Jl-v«  J  0  -ll-v  7  “ 

The  substitutions 


V  =  cos© 

( A— 12) 

V7  =  cos©7 

( A— 13) 

transform  (A— 11)  to 


U4 


ir/2 

I 

!  d©  cos© 

J0 


tt/2 

Jd©7cos©' 
0 


log | cos©-cos©’ | 


( A— 14 ) 


Substituting  C14,  eq.(A-2)3 

®  2 

log | cos©-cos© 7 |  =  -log2  -  Z  - cos(m©)  cos(m©7)  (A-15) 

m=l  m 
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S8 

J/i 

ft 


into  (A-14)  and  interchanging  the  order  of  summation  and 
integration*  we  obatin 

ir/a  tr/2 


A 

I 

U*  =  -log2  <  cos©  d©) 

J0 


which  first  reduces  to 


rr*  cd  1 

Ut  =  — log2 - 2  - 

8  m=2  2m 


®  2 

-  2  -  < 

m=  1  m  ! 


cos©  cos<m©)  d©  ) 
0 


( A-16) 


sin<  (m+1  )tt/2)  sin(  (m-1  )tt/2) 


m+1 


m-1 


( A— 17) 


and  finally  to 


Ut  =  -log2 - 2 


( A— 18 ) 


8  n=l  n(2n+l  >*(2n-l)*= 

Using  CIS,  Formulas  610.  and  610.1.1,  we  obtain 


f!  fl 

1 dv  Id v ’  1 og I v— v ’ I  = 

J  0  J0 


=  -  3/2 


< A-19) 


It  is  easy  to  show  that 


fiv  fi 

J0  J0 


V  V 

dv'  (  ■  - - 

0  Jl~vm  -J  1— v 1 


1  )  =  2 


( A— 20 ) 


Substitution  of  (A-18)-( A-20)  into  (A-3)  gives 
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Of  all  the  arguments  of  the  Hankel  functions  in  (3-89)  <■ 
only  k^C | v+v’+n-3 |  can  be  zero  at  a  point  interior  to  the 
region  of  integration.  Accordingly?  we  recast  (3-89)  as 


2knY*„  =  ( k^C)1 


0 


f1 

d 

J0 


dv-ZT-v*-  dv 9  v’  (  H0  ( k^C | v+v ’ +n-2 | > 


Ho  (k^C| v-v’+n) )  ) 


-v’+n-l  I  )  +  I, 


( A-22  > 


ri  PI 

j dv  vjdv’  v’  Ho  (k^Cjv- 

fl  v  PI  <a>  to 

dv---  ..  -  dvMHp  (  k  AC  |  v+v '  +n— 2  |  )  — H„  ( k^C | v-v ' +n | )  ) 

J  0  -i  1— v®  J  0 

ri  pi 

dv  I dv'  Ho  (k^Cjv-v’+n-l | )  +  U 

J  0  J0 


where 

1 


dv’  v’  Ho  (k^C| v+v’+n-3| ) 
0 


PI  PI 
=  dv  v  j  d 

J  0  J  0 

pi  fl 

O  =  J  dv  J  dv7  Ho  (k^Cj 


v+v '+n-3 | > 


( A— 23 ) 


( A— 24  > 


If  n  >  3j  (A-23)  and  (A-24)  are  all  right  as  they  stand 

(  S  > 

because  the  argument  of  H0  can  not  vanish  at  anv  point 


interior  to  the  region  of  integration.  If  n=2. 


ar.d  (A-24)  are  recast  as 


v  i€n  i  f- o  f 


Ia=U.+ 


!«•— U»+ 


where 


PI 


PI 

dv  v 
0  J  0 


dv'v’(H0  (kAC| v+v’-l j  >-g(k^C j v-v’-l |  ) )  (A-25) 


1 

dv 

0 


PI  <  ■ ) 

dv’  (Ho  (kACjv+v’-l j )-g(kAC| v+v'-l j ) )  (A-26) 

0 


Um  = 


1 

dv  v 
J0 


PI 

dv’v’  g(kAC| v+v’-l | ) 
0 


( A— 27 ) 
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U»  = 


‘1  PI 
dv  dv’ 

.  0  J  0 


g ( k^C | v+v’ -1 | > 


Using  CIS,  Formulas  610.,  610.1.,  610.2.,  and 

and  working  diligently,  we  obtain 


1 

dv  v 


‘1  5 

dv ’ v 1  log | v+v ’  —  1 I  = - 

.  0  16 


Substituting  (3-20)  into  (A-27)  and  using  (A-29), 
that 

1  j 

U*  =  —  +  - -  (  5  -  4 logb  ) 

4  8tt 

Substitution  of  (A-30)  into  (A-25)  gives 

fl  fl 

I3  *  dv  v  dv’vMH,  (k^Cjv+v’-l  I  )-g(k^C|  v+v’-l  |  )  1 

J  0  J  0 


1  j 

- - +  - (  5  —  4 logb)  ,  n=2 

4  Sir 

where  b  is  given  by  (A-10). 

Seeking  to  evaluate  U»  of  (A-28),  we  write 

fl  PI  3 

jdv’  dv’  log|v+v’-l|  =  - 

J0  J0  2 


To  verify  (A-32),  note  that  if  v’  is  replaced  by  l-v! 
32),  then  the  left-hand  side  of  <A-32)  reduces  to 
(A— 19).  Substituting  (3—20)  into  (A-28)  and  using  (A- 


obtain 


j 

U»  =  1  +  -  (  3  -  2 logb  ) 

TT 


Substitution  of  (A-33)  into  (A-26)  gives 


(A-28) 

610.3.1 

(A-29) 

we  find 

(A-30) 

( A— 31 ) 

( A-32 ) 

in  ( A— 
that  of 
32  > ,  we 

(A-33) 
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•1  fl  •  <«» 

=  dv  dv  * ( H„  (k<lC|  v+v’-l  |  )-g(k^C|  v+v’-l  |  )  ) 

0  0 


+  1  +  —  (  3  —  Elogb  )  ,  n=2 


( A— 34 ) 


Now*  2kqYi„  of  (3-89)  is  given  by  ( A-22 )  where  19  and  I**  are 


given  by  (A— 31)  and  (A— 34)  for  n=2  and  by  (A— 23)  and  (A— 24) 


for  n=3j4j . . . ,  M— 2. 


Of  all  the  arguments  of  the  Hankel  functions  in  (3-90), 


ily  k^Cj v+v’+M-4|  can  be  zero  at  a  point  interior  to  the 


region  of  integration.  Accordingly,  we  recast  (3-90)  as 


2krjY»  rrt—i  =  (MC)  Jdv 


1  „ _ PI 

dvJl-v*  J  d 

0  J  0 


dv'Jl-v’«=  H0  (  k^C  J  v+v ’  +M-2  J  ) 
0 


+  la  +  2 


*1  PI 
dv  v  d 

i0  t  0 


dv’4l-v’«  H„  (  kj^C  |  v+v ’  +M— 3  |  ) 
0 


'i  v  ri  v’ 

dv  -  1 dv ’ -  H„  ( k£C | v+v ’ +M— 2 | > 

0  -J  1— v®  J  0  Jl-v’  = 

PI  PI  v* 

I*  -  2  dv  dv' -  Ho  ( k^C  j  v+v '  +M-3  j  ) 

J0  J0  -il-v’“ 


( A— 35 ) 


where 


PI  PI 

[s  =  dv  v  dv’ v’  Ho 

J0  J  0 

r  r1 

to  =  dv  dv’  Ho  <  k^l 
J0  J0 


dv’ v’  Ho  ( k^C | v+v ’ +M— 4 | ) 
0 


( A— 36 ) 


C I v+v ’ +M-4 I ) 


( A-37 ) 


If  H  >  4,  (A-36)  and  (A-37)  are  all  right  as  they  stand 


because  the  argument  of  Ho  can  not  vanish  at  any  point 


KiUftUliUfWMUjtV 


.ir  I  .  •  ) 


interior  to  the  region  of  integration.  If  M=3 *  then  the 
right-hand  sides  of  (A-36)  and  (A-37)  are  the  same  as  those 
of  (A-23)  and  (A— 24)  with  n=2.  Hence?  from  (A-31)  and  (A- 
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34 ) ,  we  have 


la  — 


PI 

dv  v 


PI 


J0  J0 

1  j 

+  —  4-  - 

4  8ir 


dv’  v ’ ( Ho  (k^C| v+v’-l I )-g<k£C| v+v’-l | )  ) 


(  5  -  41ogb  )  ,M=3 


( A— 38 ) 


1 4.  = 


ri  fl 

dv 
0 


dv’  (Ho  (kiCj v+v’-l j )-g(k^C| v+v’-l | )  ) 

0 


1  +  -  (  3  -  21ogb  ) 

W 


,M=3 


( A— 39 ) 


Now,  2kqYi ,M-i  of  (3-90)  is  given  by  (A-35)  where  la  and  I6 
are  given  by  (A— 38)  and  (A— 39)  for  M=3  and  by  (A— 36)  and  (A— 
37)  for  M=4,5, - 


We  rewrite  2ki-|Ymr.  of  (3-91)  as 


\ ro=2 j  3  9 • > - « M— 2 


2kr)Ym„  =  I-r  +  la  +  I-r 


( A— 40 ) 


n-2,3, . . . ,M-2 


where 


PI  PI 
=  dv  dv’ 

J0  J  0 


(  <k4C)  vv’-  1 ) ( H0  ( k£C | v-v ’ +m-n | ) 

« m  > 

+  Ho  ( k^C | v-v ’ +n— m | )  ) 


( A— 41 ) 


la  = 


dv 


\I  YJ  O 


1  a  <  a » 

dv’(  <k^C)  vv’+  1)  Ho  ( k«2C  j  v+v  ’  +m— n— 2  j  ) 

0 

( A-42) 


mmmmmmma 


I«»  =  I  dv  |dv’(  (k^C)  vv'+  1)  H0  ( k^C  I  v+v  J+n-m-S  I  ) 

J  0  J  0 


( A— 43 ) 


If  m  k  rij  (A-41)  is  all  right  as  it  stands  because 
neither  of  the  arguments  of  its  Hankel  functions  can  vanish 
at  any  point  interior  to  the  region  of  integration.  If  m=n, 
then  (A— 41)  is  recast  as 


PI  PI 

1-,  =2U«.  +2  dv  dvMtk^C)  vv’-l)(H0  (  k/Jc  |  v-v  ’  |  )  -g  (  k^C  |  v-v '  |  )  ) 

J  0  J  0 

( A— 44 ) 


where 


’1  PI  e 

U*.  =  dv  dv’  ( (k^C)  vv’  -  1)  g ( k JC  j  v— v ’ | )  (A-45) 

J  0  J  0 


Using  (A— 7)  and  <A-19)  to  evaluate  (A-45)?  we  obtain 


(k^C)®  j2  (k^C)* 

- - 1  —  - (  <  -  —  1  )  logb 

4  it  4 


Substitution  of  (A— 46)  into  (A-44)  gives 


7  2  3 

- — (k^C)  +  -  ) 

16  2 

( A— 46 ) 


’1  PI 
=  2  dv  d 
*  0  \  0 


2  dv  dvM(kAC)  vv’-l)(H0  (  k/jC  |  v-v J  |  >  -g  (  k^C  |  v-v  *  I  )  ) 

J0  J0 

( A— 47 ) 

(k^C)®  j  s  7 

— —  -  2  -  -  (  6-  4 logb  -  (k£C)  (-  -logb)  )  , m=n 

2  TT  4 


If  m-n  41 ,  <  A— 42 )  is  all  right  as  it  stands  because  the 

argument  of  H»  can  not  vanish  at  any  point  interior  to  the 
region  of  integration.  If  m-n=l »  then  (A-42)  is  recast  as 
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i 


*1 


1 


Ia=U»+ 


where 


‘1  ri 

dv  dv ' 

0  J  0 


7((kAC)  vv’+l)(Ho  <MC|v+v’-l|  )-g(k4C|v+v'-l  |  >  ) 


Us  =  (k£C)  Ue  +  Us 


( A— 48 ) 


( A— 49 ) 


in  which  Ua  and  Ua  are  given  by  (A-S7)  and  (A-28).  In  view 
of  (A— 30)  and  (A— 33)*  substitution  of  ( A— 49 )  into  (A— 48) 
gives 


’1  ri 

—  dv  dv’ 
0  0 


((k^C)  vv'  +  lXHo,  <kAC|v+v’-l|  >-g(k4C|v+v’-lj  >) 


*  1  J  3 

+  <kAC)  (  -  + - <  5  — 41ogb  >  )  +  1  +  -  (3-2logb) 

4  Bit  n 


? m-n= 1 


( A— 50 ) 


If  n— m  1?  (A-43)  is  all  right  as  it  stands  because  the 
<  at » 

argument  of  Ho  can  not  vanish  at  any  point  interior  to  the 
region  of  integration.  If  n— m=l?  then  la  of  (A-43)  is  what 
la  of  (A-42)  would  be  if  m-n=l.  Hence? 


Ia  =  I, 


?  n-m  =1 


(A-51 ) 


where  la  is  given  by  (A-50). 


Now?  2krjYm„  of  (3-91)  is  given  by  (A-40)  where  la  is 
given  by  (A-41)  for  m  ^n  and  by  (A-47)  for  m=n?  Ia  is  given 
by  (A-42)  for  m-n  ^  1  and  by  (A-50)  for  m— n=l?  and  la  is 
given  by  (A-43)  for  n-m  ^  1  and  by  (A-51)  for  n-m=l . 
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APPENDIX  B 


In  this  appendix,  we  summarize  the  scattering  method,  a 
method  in  which  the  aperture  remains  open  and  only  electric 
current  is  involved  in  the  solution.  In  the  original  problem 
in  Fig.l,  the  impressed  sources  radiate  in  the  presence  o-f 
the  coducting  bodies.  The  resulting  -Field,  denoted  by 
can  be  viewed  as  the  sum  o-f  the  -Field, 
denoted  by  Eloe,  which  would  exist  i-f  the  impressed  sources 
radiated  in  -Free  space  and  the  field,  denoted  by  E(J“,Q), 
produced  by  the  eletric  current  J*  induced  on  the  perfect 
conductors  radiating  in  free  space. 

The  boundary  condition  is  that 


A  n 


E* 

+  E 

Q)  = 

:  Q 

on  S«,  , 

and 

(B— 1  > 

S, 

where  t 

denotes 

the 

tangential 

componet  on 

the 

conducting 

surface  and 

S-, 

s*. 

and 

5«  are 

the  surfaces 

of  the  perfect 

conductor 

5- 

(See 

Fig. 

1) 

(B-l) 

states  that 

the 

tangential 

componet 

of 

the 

el ectri c 

field 

vanishes  on 

the 

surface  of 

the  perfect  conductors.  Rearranging  (B-l),  we  obtain 

1  nc 

-  E* (J-,0)  =  E«  (B— 2) 

To  emphasize,  we  repeat  that  Etr,*=  is  the  electric  field  that 


would  exist  if  the  impressed  sources  radiated  in  free  space, 
and  E(J*,0)  is  the  electric  field  due  to  J*  radiating  in 
free  space  as  well.  E(J*,0)  is  usually  called  the  scattered 


field.  Thus,  we  refer  to  the  method  developed  here  as  the 
scattering  method. 

In  the  following,  we  obtain  J*  by  the  moment  method 
C51.  Let  J*  be  approximated  by 

N  . 

J-  =  I  I„  Jj  (B— 3) 

j  =  l 

where  i J j }  are  expansion  functions  which  must  be  chosen 
tangent  to  the  conducting  surfaces,  flj>  are  coefficients, 
and  N  is  the  number  of  expansion  functions  used  to 
approximate  J“. 

Substituting  (B-3)  into  (B— 2) ,  we  obtain 

N  * 

~  £  Ij  =  E«;  (B-4) 

J  =  1 

m 

where  E(Jj,0)  is  the  electric  field  due  to  the  j tri  expansion 
funciton  J_,  radiating  in  free  space.  Taking  the  summetric 

m 

product  of  (B-4)  with  J x ,  i=l,2,...,N,  we  obtain 

CZ-:  I"  =  V  (B-5) 


where  the  ij**’  element  of  CZ*3  is  given  by  (2-24)  with  the 

m 

superscript  a  replaced  by  s,  the  element  of  I»  is  Ij, 

*  -> 

and  the  i**1  element  Vi  of  V*  is  given  by 


•  •  1  nc 

Vi  =  ^  Ji ,  E*  > 


(B-6) 


After  1“  is  obtained  by  solving  the  set  of  linear  equations 


in  (B-5),  the  fields  can  be  written  as 
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£ 


*1*1 


'jfj 


yl 

*.ii 


.  N  . 

H  ( J 1  "*«* , M1  "•*» )  =  H‘"«  +  Z  Ij  H(Jj,0) 

j  =  l 


(B— 8) 


where  each  field  quantity  on  the  right-hand  sides  of  (B-7) 
and  (B-8)  is  due  to  its  source  radiating  in  free  space. 


The  formulation  here  is  indeed  simple.  Unfortunately, 
it  will  fail  to  accurately  determine  the  field  inside  the 
cavity  formed  by  the  walls  of  a  conducting  body  if  this 
cavity  has  a  very  small  aperture.  This  is  due  to  the  fact 
that  the  field  inside  the  cavity,  being  the  sum  of  the 
incident  field  and  the  field  due  to  J*,  will  be  very  small. 
The  field  due  to  J“  nearly  cancels  the  incident  field. 
Hence,  a  small  percentage  error  in  J*  will  give  rise  to  a 
large  percentage  error  in  the  field  inside  the  cavity. 


•vVvSVm 
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